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CHAPTER I: EXECUTIVE SUMMARY 


Section 1: Need for Study 

Various Investigations and model calculations yield evidence that 
thermal Infrared remote sensors are useful In assessing water-related 
surface and near-surface land conditions. Soil moisture Is highly 
dynamic having extreme spatial and temporal variability. This charac- 
teristic makes ground point- sampling observation programs expensive and 
relatively Ineffective for the sequential and synoptici yet detailed, 
Information required In management trograms. Remote sensing systems have 
the capability to provide these data characteristics If the Information 
content of the data can be quantitatively used to assess moisture related 
ground variables. 

The agriculturalist has an Interest In soil moisture which Is 
available to the growing plants. This varies with season and crop but 
ranges from surface moisture to profile soil moisture at depths of a 
meter or more. Conversely, the hydrologist may have an Interest only 
In the surface few centimeters of the land surface. Of special Interest 
In South Dakota Is the depth to the zone of saturation. Where this 
depth Is within a few centimeters to approximately one meter, the soils 
may have a limited rooting depth and may tend to develop soil salinity 
caused by upward mobility of the water table moisture with subsequent 
evaporation and deposition of salts. The contamination of water table 
water associated with iw« Ing of the soil horizons has a higher poten- 
tial If the water table Is near the land surface. 


The depth to the saturated zone provides an Indication of the water 
budget of surficlal aquifers. Where discharge through evapotranspi ration 
of phreatrophytes or of Irrigated crops and the domestic consumption of 
water In combination with other loss factors exceeds recharge* water 
tables will lower. The Importance of understanding aquifer water budgets 
Is acknowledged since tedious ground observation programs are already In 
place In eastern South Dakota. A remote sensing program could potentially 
aid these programs for appropriate placement of sounding wells or Infer- 
ring contours among the point observations for final mapping. 

Ter many areas In Sojth Dakota and other regions* monitoring the 
depth to the saturated zone* where It occurs In the top 2-3 meters of 
the land surface* would provide a valuable Input into water and land 
management. Remote sensing approaches to a moi'ltcring program would 
use surface or near-surface measurements (I.e. thermal or microwave) to 
Infer certain water-related properties. Therefore, changes in the energy 
budget of the surface which result from a variety of factors must be 
evaluated and understood before thermal -Infrared remote sensing data can 
provide the desired Information. The shading effects of the crop canopy, 
water use by evapotranspi rati on* incomplete canopy cover, temperature 
indicators associated with variations of thermal diffusivity* slope, etc., 
are all candidate topics which must be addressed and understood to fully 
utilize the data derived from thermography. 

The objectives pursued under this Investigation were to: 

(1) Develop a finite-difference simulation model for specific 
application to HCMM data which simulates and predicts the 
thermal regime associated with occurrence of shallow groundwater. 


(2) Evaluate HCMM data for locating regions of shallow groundwater 
and estimating the depth to the water table. 

Section 2: Report Format 

Detailed results are presented as various chapters. Each chapter 
Includes published or submitted paoors which Identify results relevant 
to the dual objectives. An additional topic beyond the original objec- 
tives was Investigated which was to determine the uses of HCMM or 
similar data In soil geography studies. 

The specific background » procedures* and discussion are presented 
within each chapter. A generalized summary of the Investigation results 
Is presented In the following section. Since potential applications of 
HCMM-type data have many questions to be resolved* any number of valid 
experimental topics could be addressed. Many were pursued under this 
program using ground or ali^craft observatlor.<; since It was questionable 
during the Investigation If appropriate HCMM data would be received in a 
timely manner to use the data for d.:alys1s. Therefore, results and pub- 
lications were fragmented ind are presented In entirety in this report 
since procedures, test sites, etc. were not consistent for all investi- 
gations. 

Section 3: Suimary of Significant Results 

Section 3.1. - Introduction 

Landforms associated with areas of shallow groundwater In glaciated 
eastern South Dakota can be recognized using standard image Interpretation 
techniques of multidate Landsat data. Therefore, this Investigation was 
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to determine, once these regions were Identified and mar)ed. If synoptic 
thermal data from HCMM could be used to assess the depth to the water 
table. Since the thermal signal measured by HCMM is a function of 
emitted radiation from the atmosphere and from the land surface, both 
must be considered when relating HCW1 data to land or water features. 

The approach to reduce atmospheric variation used spatial differences 
within dates for analysis. No further attempts were utilized to reduce 
the Influence of the atmosphere on the thermal signal even though the 
authors realize this as a significant problem in quantitative rerote 
sensing using thermography. Therefore, most equations and statistical 
relationships presented In this rtport require that some form of 
quantitative field data are required to Implement 3ny of the algorithms. 

Since the composite "surface" radiance is that which Is sensed 
from the remote sensor, evaluations of the effects of vegetation cover 
^nd »oil moisture were investigated. Sections 1-5 of Chapter II detail 
each specific stud; . In general, the topics 1r Ijded the use of thermal 
data for assessing actual soil temperatures under conditions of Incom- 
plete and complete canopy covers and the assessment of near-surface soil 
moisture. Chapter III follows with the relationship of ground, aircraft, 
and HCMM data to the depth- to-wa ter tables. Since tiie program had 
considerable delays in dissemination of HCMM data, another objective was 
pursued to evaluate HCMM for use in operational soil surveys. Chapter IV 
clearly Illustrates and demonstrates examples where HCMM or similar data 
could provide input to the field soil surveyor. The final Chapter V 
defines a model and its results for developing a theoretical base to 


5 


understand thermal Infrared potentialities In water resources monitoring. 

Section 3.2 • Qualitative Monitoring of Wet Dry Surfaces 

In Its simplest form, daytime thermography should provide a qualita- 
tive mapping tool to locate regions which have a water saturated surface 
versus a surface where the soil moisture Is at less than saturation. 

The effects of higher thermal Inertia and 1ncreai';;d partitioning of in- 
oiimlng energy Into latent heat (Increased evapotranspi ration^ result In 
cooler surface temperatures. Chapter II, Section 1 (1 1-1) and Figures 
1-3 Illustrate the phenomena. The predawn Image In Figure 1 demonstrates 
that the area of a cool anomaly presented ^n Figure 2 was probably both 
a function of differences In evapotranspiratlon and of thermal Inertia. 
The warmer apparj ■ tempe*'"tures predawn results from a thermal Inertia 
change of Increased moisture when compared to surrounding regions. The 
cooler apparent temperatures of the saturated area during daytime hours 
results from both thermal Inertia variations and Increased evapotrans- 
piratlon. 

Section 3.3 - Thermography of Incomplete Canopy Cover 

The remotely measured land surface emission Includes radiation 
emitted from both soil and vegetated surfaces. If the Intended applica- 
tion requires either soil or vegetation temperaturei Indepenlently, a 
method to separate and calibrate the contribution of each emitting com- 
ponent must be developed. For the water table monitoring objective, 
the need Is to measure soil surface temperature. If the application Is 
to monitor transpiration of plants as an Indicator of moisture or other 
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Figure 1. Photographic enlargement of a May 14, 1978 night thermal infra 
red image (scene I.D. A-AOOl 8-08420) showing a high soil mois- 
ture area (arrows) in southeastern South Dakota. Thermal 
Inertia is apparently causing a warmer signature during the 
predawn cooling portion of the diurnal cycle. Dark is cool. 
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Figure 2 . Photographic enlargement of a May 15, 1978 day thermal Infrared 
Image (scene I.D. A-A0029-19575) showing a high soil moisture 
area (arrows) In southeastern South Dakota. Tl,e area Is 
approximately that of Figures 1 and 3. Increased evapotranspira 
tion and Increased thermal Inertia are probably responsible 
for the cooler apparent temperatures of the wetter area. Dark 
Is cool. 
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Figure 3. Photographic enlargement of a May 13, 1978 Landsat MSS7 
image (scene I.D. E-21207-16083) of the same approximate 
area as Figures 1 and 2. The area i>^ question does not 
fiave standing surface water but has ^ther bare or vege- 
tated fields as the emitting surface. 




stre<«s, the vegetation component must be estimated. Section 2 in Chapter 
II presents an algorithm to extract canopy temperature from radiometric 
measurements at Incomplete cover. 

A ground-based study was conducted on a developing barley canopy 
In which canopy temperatures were measured “Hh an Infrared radiometer 
at 30° from horizontal. Contact soil temperatures at a 1-nm depth were 
measured using thermocouples. Composite radiometric temprratures which 
Included radiance contributions both from the soil and from the barley 
vegetation were measured with a radiation thermometer at 2-m above the 
canopy. This measurement simulates that of vertical remote sensors. 
Emissivlty and sky Irradlance corrections were applied to the data. The 
measut^ments were acquired at approximately 1330 1ST ( -.'!M overpass time) 
throughout the barley growing season. The canopy covers seasonably 
varied from 30^!o to 90% (*LAI 0.3 to 3.2). 

The radiometric temperatures of the composite radiating surface, 
Including both soil and crop canopy, were 0.5 to 11. F C higher than 
canopy temperatures alone. Surface soil temperatures were 1.5 to 20.0 C 
higher than canopy temperatures. Correlation between measured canopy 
temperatures and composite temperatures was not statistically signifi- 
cant. 

The difference between composite and canopy temperatures was highly 
correlated with percent canopy cover (r ■ 0.52). Therefore, an approxi- 
mation equation was applied which separated the radiance contributions 
of the bare soil and the canopy. The equation incorporated terms for 
canopy and soil emissivltles, percent crop cover, longwave sky irradlance, 
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and composite temperatures. Applying the correction resulted In a 
2 

r • 0.88 for relating predicted versus observed values for estimating 
the temperature of the canopy. Differences between nr«i<jictsd and ob- 
served values ranged from -1.84 to 2.5 C for the data acquired throughout 
the barley growing season. 

The significance of the procedure was. even at low canopy covers 
estimates of canopy or soil temperatures can be successfully measured 
using remote sensors. Figure 4 summarize} the data for the barley study. 
Remote sensing estimates of percent cover and composite temperatures 
are available with existing technology. When emissivlty corrections 
were not applied, prediction accuracy varied >t1th percent cover. The 
largest errors occurred at low cai.opy cover. Canopy emissivlty was 
measured as 0.98 and was assumed constant throughout the growing season. 
Soil emissivlty varied with water content. Measured values ranged be- 
tween 0.95 and 0.97. Wider ranges of soil emissivitles can be experi- 
enced since the value Is soil specific. If emissivitles were assumed 
to be one for the calculations In this field experiment, predicted 
canopy temperatures would range from 6.4 C lower to 1.7 higher than 
actual observed values. Where measured emissivitles were used In the 
equation but the longwave sky Irradlance was Ignored, differences be- 
tween predicted and observed canopy temperatures ranged from 0.<3 to 
10.7 C. 

The significance of the study was that canopy temperatures at In- 
complete covers can be estimated using a nadlr-viawlr.g radiometer If 
appropriate considerations are given to soil background radiance, 


n 



Figure 4. Comparison of predicted and observed canopy temperatures 
where the predicted temperatures utilized terms of canopy 
and soil emissivities, percent crop cover, longwave sky 
irradiance. and composite radiometric temperatures (or 
remotely sensed temperatures). 
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emissivltles. and sky Irradlanre. Techniques are available for measur- 
ing or estimating some of the inputs required for the canopy model. 

Various models which have been developed for simulating bare surface 
soil temperatures can feasibly be extended to soil surfaces including 
crop canopies. 

Section 3.4 - Diurnal Surface Temperature Change as a Measure of Near- 
Surface Soil Moisture 

The basic assumption in using thermography to estimate near-surface 
soil moisture is that surface temperatures or their changes in time are 
related to the thermal properties of soils in a predictable fashion when 
other components of the energy budget are not accounted for or can be 
accounted by use of auxiliary data which are easily available. The 
paper presented in Chapter II, Section 3 describes in detail the experi- 
mental framework to determine if aT of the surface relates to near-surface 
nioisturfi and if this relationship is maintained with the variability 
introduced throughout a barley growing season having variations in canopy 
cover. The same rainfed barley field as described in the previous sec- 
tion was instrumented and measured on 22 dates over a 45-day period. 

Volumetric soil water content of the 0-4 cm layer was the variable 
to be predicted by thermography. The field data were cfillected through- 
out a barley growing season where percent cover ranged from 30% to 90% 
and leaf area index from 0.3 to 3.2. Diurnal changes in 1-mm soil 
temperatures between 1330 and 0230 LST (aT^) were used in an exponential 
^ equation to predict soil moisture content of the 0-4 cm soil layer. An 

^ - j) 

r = 0.81 was derived for the equation which incorporated data from the 
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Figure 5. Relationship of the difference (Ms) between soil surface 
temperatures measured at 1330 and 0230 LSI and the average 
24 hr volumetric soil water content in the 0 to 4 cm layer 
of the profile. 


total season. The resulting equation and plot Is Illustrated In Figure 
5. Attempts to normalize among dates using air temperatures provided 
no significant Improvement to the results. 

The quantitative relationship Is significant In Illustrating that 
If an estimate of actual dally surface temperature difference Is available. 
It Is statistically highly significant In direct relation to near-surface 
moisture. The diversity In canopy cover, a1i 'temperature and humidity, 
and other energy budget terms did not mask the relationship. 

Section 3.5 - Surface Soil Temperature and Moisture Estimates Using 
Radlometry. 

The same barley canopy and experimental data were evaluated to 
determine If radiometric temperatures of "composite” land surfaces (In- 
cluding both canopy and bare soil) could be used to quantitatively re- 
late to the 0-4 cm soil moisture. The actual 0230 LST composite temp- 
eratures ranged from 1.1 - 2.2 C higher than canopy temperatures (ex- 
cluding soil backg**ound) . This predawn background soil Irradlance 
was typically higher than the canopy temperature. The surface soil 
temperatures were higher by 1.1 - 5.4 C than canopy temperatures. This 
trend remains during the 1330 LST measurements. The composite radio- 
metric temperatures were 0.5 - 17.0 C higher and surface soil tempera- 
tures 1.5 - 20.0 C higher than temperatures of the vegetation in the 
canopy. Therefore, changing canopy cover affects the ability to esti- 
mate the actual soil surface temperature which Is the temperature that 
is related to near-surface moisture. 
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Regression equations were developed using the seasonal data of the 
barley canopy to relate composite (or remotely sensed temperatures) to 
the actual soil surface temperatures. For the 0230 LSI measurement time* 
the resulting equation Included modifying the thermography measurement 
with National Weather Service dally minimum air temperature. For the 
1330 LST measurement, the equation Included a term modifying the 
thermography with an exponential term of fraction of vegetation cover. 

The two equations had R *0.78 and r*" • 0.86, respectively. "Igure 6 
Illustrates the corrected data after application of the canopy cover 
correction equations. 

The equations were derived empirically and should be questioned and 
evaluated over a larger set of environmental variables than were Included 
In this original data set. However, the need for some type of correction 
of the thermal radiance measure Is required to relate to actual soil 
surface conditions where partial or Incomplete cover by a crop canopy 
Is present. 

The equations presented In this study Illustrate that easily derived 
field properties such as air temperature or percent canopy cover may pro- 
vide the necessary Information. The equations developed were limited 
In geographic extent but did encompass data throughout a complete 
growing cycle of one crop. 

Section 3.6 - Aerial Thermography as an Estimate of Near-surface Soil 
Moisture Under a Variety of Land Covers. 

The correction equations developed from the barley study were 
J applied to an aircraft study of a larger region In eastern South Dakota 

which Included a diversity of canopy covers (50%-95% cover Including 
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Figure 6. Comparison of predicted and observed values of surface soil 
temperature at 0230 LST and 1330 LST. 
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corn, soybeans, millet, and pasture) and differences In soil textures 
and other soil properties (Chapter II, Section 3). The day/night air- 
craft thermal data were collected on only one date. Soil nolsture was 
estimated as percent of field capacity to reduce variations associated 
with soils. The average difference between measured and predicted 
values was 1.6% of field capacity. Figure 7 provides a plot of the 
measured versus predicted soil water contents. 

The significance of the study was that the correction equations 
were applicable to a wide range of crop cover conditions. This range 
Is normally prevalent when applying the techniques over large agricul- 
tural areas where soil moisture monitoring programs might be established. 

Section 3.7 - HCMM ThenBography as an Estimate of Near-surface 
Moisture Under a Variety of Land Covers 

Satellite applications have the advantage of covering large regions 
in a short time period. When the theme to be estimated Is confused with 
other scene variables, the complexity caused by heterogeneity of the 
terrain normally Increases as larger areas are ooserved. The same 
correction procedures were tested using data as thermography covering 
multiple dates over a sizable area In eastern South Dakota having a com- 
plex agricultural landscape. 

HCMM data were treated with correction equations similar to those 
discussed In the previous paragraph (Chapter II, Section 5). Soil 
moisture data from the 0-4 cm soil layer were acquired from 23 fields 
having a variety of crop covers for four dates ranging from June 
through early September. The HCMM radiometric temperatures without 
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Figure 7. Comparisc.. of predicted and observed values of 24 h average 

soil water content In the 0 to 4 cm layer of the soil profile. 
Aircraft thermal scanner data were used In prediction after 
applying correction equations for canopy variations. 
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correction for land cover did not significantly relate to soil moisture. 
However, when actual surface soil temperatures were empirically estimated 
from HCJ#I radiometric temperatures using land cover corrections, a highly 
significant correlation (r ■ 0.74) was obtained betw4>en the estimated 
soil temperatures and 0-4 cm soil moisture content. Figure 8 Illustrates 
the derived relationship. 

Sect. on 3.8 - Depth to Water Table 

If the water table or zone of saturation Is at a greater depth than 
the diurnal damping depth, the amplitude of the diurnal temperature curve 
1: not affected but a shift In the curve up or down In absolute magnitude 
will be noticed where the depth to saturation Is within the depth of 
annual temperature variations. Many Investigations, as reported In 
Chapter HI. Sections 1 and 2. have demonstrated that 50-cm deep soil 
temperatures will vary with depth to water table. The saturated zone 
acts as a heat sink or source which delays warming Id the summer and 
retards cooling In the winter. 

The Big Sioux River Basin In eastern South Dakota Is a surficlal 
shallow water table aquifer of glacial origin. The land use Is predomi- 
nantly agriculture including small grains, pasture, hayland. and row 
crops. A considerable difference In the thermography can be expected 
as associated with land cover. Contact temperature measurements at a 
50-cm depth were acquired to determine If a significant relationship 
between the soil profile temperatures and depth to water table was 
retained even under varying land cover (Figure 9). Similar to 
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Figure 8. Relationship between predicted soil surface temperatures (T ) 
ard soil water content In the 0-4 cm laye-s of the soil 
profile. Field soil water contents are an average of three 
measurements. 
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Figure 9. Relationship of 50-cm soil temperatures to depth to ground- 
water for row crops and pasture. Temperatjres were measured 
in the Big Sioux River Basin during daylight hours on 
September 5-7, 1978. 
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observations •"eported In the literature, a significant relationship be- 
tween measurea JO-cm soil temperatures and depth to water table was 
observed for water tables in the Big Sioux Basin of eastern South 
Dakota. Three meters was the depth limitation of a significant correlation. 
These correlations were observed during early September near the annual 
maximum downward temperature gradient. The saturated soil materials 
act as a heat sink as the annual temperature mark progresses into the 
summer in the northern latitude of Sou\h Dakota. This creates cooler 
land surfaces over the water saturated zone. However, both soil moisture 
and vegetation cover may obscure any surface expression of this tempera- 
ature anomaly. 

The range of 50-cm temperatures was about 5 C for the data acquired 
under this investigation. The 50-cm temperature anomaly predicted by 
modeling and f*eH observations has previously been shown in the litera- 
ture to be 1-2 C as associated with water tables. However, data points 
in this study were in soil profiles having a variety of crop cover 
conditions on the surface which also affected the 50-cm temperatures. 

High soil moisture and shallow groundwater tend to affect surface tempera- 
tures in the same direction - the surface remains cooler during the day 
{in periods of maximum downward temperature gradients - August or 
September) - but in opposite directions at night. At night, the high 
soil moisture areas should remain warmer when considering thermal inertia, 
and the shallow groundwater areas should remain cooler. 

Previous investigations have shown that aircraft radiometric 
temperatures were related to water table properties. Predawn HCMM 
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Imagery was effectively used with cooler anomalies delineating the 
alluvial basin when the water tables were shallow In Figure 10. Land* 
form recognition from Landsat or aircraft Image data was also used to 
locate these potential areas. During the day, the HCMM thermal data 
did not reveal the locations of the basin as Illustrated In Figure 11. 

The landforms (alluvium and terraces) where high water tables are 
commonly located were delineated using a combination of Landsat and 
HCM1 predawn data. Since the depth to the water table had significant 
correlations to the 50-cm temperatures, HCMM data were evaluated to de- 
termine If the radiometric surface temperatures were correlated with 
water table depths. Uncorrected HCMM data were evaluated for estimating 
water table depths within the basin for four dates from June through 
early September. No significant relationship was found. The data 
were corrected using the empirical vegetation transformation developed 
under the previous soil moisture studies to reduce the effect of the 
crop canopies when estimating surface temperatures. Each date was 
analyzed separately. The correlations Improved from 0.59 to 0.8 as 
the season progressed from June 5 through September 4. The September 4 
date Is plotted In Figure 12 which relates corrected HCMM data to water 
table depth. 

This parallels model predictions. The maximum anomaly should occur 

In the season of greatest downward temperature gradient. Since both 

‘.oil moisture and water table depth can affect surface soil temperature. 

a multiple regression analysis was conducted to determine the proportion 

of the HCWi radiometric temperature which could be accounted for with 

2 

these two variables. The resulting R * 0.87 was for the September 4 
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Figure 10, Photographic enlargement of an August 29, 1978 night HCMM 
thennal infrared image (scene ID A-A-1 25-08340) showing 
the Big Sioux Basin. Note that the Basin appears cooler 
than the surrounding areas, due primarily to the heat sink 
produced by shallow groundwater within the Basin. 
(Approximate scale 1:1,000,000; dark is cool) 
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Fioure 11- A September 4, 1978 photographic enlargement of a HCMM day 
^ ^ thermal Infrared image (scene ID A-A0131-19420) of the same 

jl area shown in Figure 10. Note that the Big Sioux River Basin 

is not visible because of emittance variation ?iSSOciated with 
land use (Approximate scale 1:1,000,000; dark is cool). 
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Figure 12. Relationship between predicted soil surface temperature 
and water table depth on September 4, 1978. 
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2 

date. This compared to r •0.64 when relating only water table to HCMM 
thermography. Please note that depth to the water often Influences soil 
moisture. 

Data such as HCMM can be beneficial In guiding well drilling pro- 
grams for monitoring and mapping water table depths where the water 
tables range from 0 to 3 m. Land cover masks the thermal anomaly 
associated with water tables but the effect can be reduced by using pre- 
dawn data or by accounting for land cover In the prediction equation. 
Since the thermography signal Is a function of many variables Including 
the heat sink (or source), land cover, and soil moisture, results of 
the Investigation revealed that If the latter two are Introduced Into 
equations to assess depth to water table. Improved correlations can be 
der^ ved. 

Section 3.9 - Soil Geography 

Soil moisture, vegetation cover, elevation, slope aspect, soil 
texture, and many additional variables affect land surface temperatures. 
Many of these variables are useful to the soil surveyor who Is classi- 
fying soils and estimating land capabilities. Image interpretation of 
HCMM thermal data of South Dakota was conducted to determine Its use 
for the soil surveyor. Results were that elevated areas, soil textural 
differences, and slope aspect were apparent in the imagery. Warmer 
radiometric predawn temperatures of north-facing slopes were evident. 

The data provide a pictoral Illustration of how soil moisture and vege- 
tation variations are associated with aspect of slope. North-facing 
slope aspect reduces solar Insolation to the surface resulting in less 


^wapotranspl ration and higher soil moisture. HCMM or similar data can 
provide this Information by estimating the net effect of slope aspect 
as It relates to vegetation growth and production (Chapter IV). Simi- 
lar observations were made for elevation differences between areas where 
parent materials and soils appear similar but noted differences In soil 
productivity are present. Figures 13, 14 and 15 Illustrate the principles. 

Section 3.10 - Model Development 

Since the surface temperatures as measured radlometrically by remote 
sensors can vary with many factors, a finite-difference simulation model 
was developed to Isolate variables and estimate their magnitude of effect. 
The concept was to measure profile temperatures. Insolation, etc. at one 
site and correlate the differences In variables which are present at a 
second site and are causing the differences In radiometric temperature. 
Using this spatial differencing approach also reduces the effects of 
atmospheric disturbances. The model (Chapter V) showed that the 50-cm 
soil temperature difference associated with water table depth differences 
between two sites was apparent as temperature differences at the surface. 
These surface differences were nearly constant 1n magnitude throughout 
the diurnal cycle. However, the difference was destroyed during pre- 
dawn In areas having dense canopy covers. Therefore, remote measurements 
can be day or night, whenever the actual soil surface temperature related 
best to the radiometric temperature as measured with remote sensors if 
dense crop cover Is not present. Soils with different moisture profiles 
differed In surface temperature variations during the diurnal cycle. 

This functional form of differences was changed less by variations in 
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Figure 13. Predawn thermal infrared (HCMM) image for an area in 
east-central South Dakota (dark is cool). 
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cinopy covtr than wit tht functional fom of tho Individual titti. 
MuUlpla pats satalllta data should b« acMuIrtd during daylight hours 
according to tha modal pradlctlons for attatting tamparatura diffarancat 
ralatad to soil molstura. Modal calculations using day minus night 
tamparaturas wara not utaful for atsatting watar tabla dapths sinca tha 
tharmal Inartlal changa was not within tha diurnal damping dapth. 

Elthar day or night data could ba utad for this ralatlonshlp. If danta 
canopy covar was prasant, tha timing of data coll action was 11ml tad to 
day pariods. 
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BACKGROUND 

Thtrmal Infrartd datuctlon and quanta M cation of naar>surfaca soil 
water content are based on relationships between surface soil tempera- 
ture and soil moisture. Diurnal variations of surfjice soil temperature 
are related to soil thermal properties and meteorological factors such 
as solar radiation, air temperature, relative humidity, wind, etc. The 

meteorological factors represent the driving force for diumai soil 

-2 -I* -1 

temperature variations. Thermal Inertia (Jm sec ^ K ), defined as 
(aC)** where x(wm”^ K'b Is thermal conductivity and C(Jm”^ K“^) is 
volumetric heat capacity, represents a sol Vs resistance to the driving 
force. Since x and C Increase with an Increase of soil moisture, the 
resulting range of surface soil temperature will decrease. 

When the soil surface Is wet, evaporation is a major factor 
controlling surface heat loss. As the surface layer dries and the soil 
water supoly cannot meet the evaporative demand, soil temperature Is 
largely Influenced by thermal Inertia. Thus, the diurnal ange of 
surface soil temperature can be an indication of soil water content. 
Idso et al. (1975) found a significant relationship between the diurnal 
range of s rface soil temperature (bare soil) and surface soil water 
content, and reported that the relationship was a function of soil 
type. Pratt and Ellyett (1979) presented a method for estimating soil 
thermal properties for changes in composition, porosity, and moisture 
content. Although temperature versus water content relationships are 
complicated by vegetation, Heilman et al. (1978) demonstrated the 
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potential for estimating near-surface soil moisture from remote 
temperature measurements of crop canopies at Incomplete cover. 

One objective of NASA's Heat Capacity Mapping Mission (HCMM) 

1' to evaluate the feasibility of using HCMM data to assess soil mois- 
ture effects by observing temperatures near the maximum and minimum of 
the diurnal temperature cycle. The satellltet which carries a two- 
channel radiometer (0.5 to 1.1 and 10.5 to 12.5 vm). collects data at 
1:30 p.m. and 2:30 a.m. local time at mid latitudes with a repeat cycle 
of 5 or 16 days depending on latitude. Spatial resolutions are 0.5 x 
0.5 km at nauir for the visible channel and 0.6 x 0.6 km at nadir for 
the thermal Infrared channel. An example of HCMM detection of a 
region of high soil moisture Is presented In the following discussion. 

DISCUSSION 

In early April 1978 heavy runoff from snownelt and Ice blockage 
caused significant flooding of alluvial areas In a portion of the Big 
Sioux River Basin In southeastern South Dakota (Figure 1). By mid-May* 
flood waters had receded, but an area of high soil moisture (at or near 
field capacity) remained. Soil moisture In the surrounding terrace 
soils was generally less than in the flood plain. 

The high moisture area appeared wanner than surrounding areas on 
May 14 HCMM night themtl imagery (Figure 2) and cooler than surroinding 
areas on May 15 HC^W day thermal Imagery (Figure 3). The temperature 
differences between alluvial and surrounding areas were probably the 
result of thermal inertia and evaporaticn differences associated with 
soil moisture differences. The high moisture area was not visible on 
Landsat Imagery (Figure 4). 
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F1g. 1. Landform map of Brookings County, South Dakota, 

showing location of alluvial soils (bottomland) of 
the Big Sioux River Basin which were flooded In 
early April 1978. 
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2. Photographic enlargement of a May 14, 1978, night thermal 
infrared image (scene ID A-A0018-08420) showing a high 
soil moisture area (arrows) in southeastern South Dakota. 
Dark is cool . 
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Fig. 3 Photoyt a^;’iic enlargement of a May 15, 1978, day thermal 
infrared image (scene ID A-A0029-1 95''5) showing a high 
soil moisture area (arrows) in southeastern South Dakota. 
Dark is cool . 
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Fig. 4. Photographic enlargement of a May 13, 1978, Landsat MSS '/ 
image (scene ID E-21207-16083) of the same area shown in 
Fig. 2 and 3. 
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Although digital data were not available at the time of the 
writing of this article to quantify radiometric temperature differences 
associated with the soil moisture differences, results presented here 
demonstrated the superiority of HCMM thermal data acquired at the appro- 
priate periods of the diurnal temperature cycle over Landsat data for 
assessing soil moisture differences. Final results from HCMM soil 
moisture Investigations currently In progress will fully evaluate the 
utility of using HCMM and similar data for evaluating soil moisture from 


space. 
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ABSTRACT 


A field study was conducted In a barley {Hordeum vutgax^e L.) canopy 
to assess the potential for extracting canopy temperature Information from 
nadir radiometric measurements at Incomplete cover. Composite tempera- 
tures consisting of emitted and reflected longwave radiation from the 
barley and the soil background were measured by a nadir-viewing Infrared 
radiometer. Canopy temperatures were measured by an Infrared radiometer 
at a 30° angle from the horizontal. Soil temperatures were measured 
with thermocouples. 

Composite temperatures were 0.5 to 11.5 C higher than canopy tempera- 
tures with the largest difference occurring at low canopy cover. The 
correlation between composite and canopy temperature for data acquired 
throughout the growing season was not significant. A model which con- 
sidered emitted radiation from both the canopy and the soil background, 
and which Included reflected longwave sky Irradlance was used to predict 
crop temperatures from nadir measurements. Predicted temperatures agreed 
with observed values (r^ ■ 0.88), and the prediction accuracy was Inde- 
pendent of canopy cover. When emissivity corrections were not applied, 
prediction accuracy varied with percent cover with largest errors occur- 
ring at low cover. Prediction accuracy also varied with canopy cover 
when appropriate emissivltles were used but sky Irradlance was ignored. 
Results Indicate that canopy temperatures can be estimated from nadir 
measurements at Incomplete cover If percent cover, soil temperature, 
soil and canopy emissivltles, and sky irradlance are known. 
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INTROnUCTION 

Remotely-sensed surface temperatures can be useful for many agricul- 
tural appllcotlons Including evapotranspiratlon modeling, soil moisture 
detection, plant stress detection, yield prediction, and Irrigation 
scheduling. Most studies which have used remote measurements have been 
restricted to bare soils or fully developed crop canopies because of the 
complexities Involved In Interpreting thermal data at less than full 
cover. 

Much of the complexity results because the remote sensing Instrument 
measures emitted and reflected radiation from vegetation and soil which 
generally have different temperatures and emissivltles. Hatfield (1979) 
reported that differences betv^en angular and vertical Infrared ther- 
mometer n»asurements of canopy temperatures were greatest at 20 to 50% 
cover and decreased as canopy density Increased. He speculated that 
differences were enhanced by emisslvity variations. Millard et al. 

(1980) found that for canopies covering at least 85% of the soil surface, 
airborne measurements of plant temperatures differed from ground measure- 
ments by less than 2 C. At 50% cover, differences were as large as 9 C. 
Investigators have shown that even at full cover thermal radiance from 
the soil surface can affect remote temperature measurements of crop 
canopies (Blad and Rosenberg, 1976). 

Incomplete plant canopies are important remote sensing targets be- 
cause of the potential benefits arising from early assessment of crop 
condition. Jackson et al. (1979) presented a model for extracting 
canopy temperature information from a composite of soil and plant 
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temperatures measured by a sensor scanning perpendicular to canopy rows. 
They found that If a critical scan angle (determined from reflectance 
measurements) was exceeded, the ten^rature obtained from the scanner 
was that of sunlit vegetation. They also found that the extraction 
process was difficult for canopies having low percent cover. 

We evaluated relationships among percent cover, soil temperature, 
and radiometric measurements of canopy temperature, and used a model of 
emitted canopy and soil background radiance and reflected sky Irradlance 
to assess the potential for extracting canopy temperatures from nadir 
radiometric measurements. Effects of neglecting emissivlty variations 
and sky Irradlance In the model calculations were evaluated. 

MATERIALS AND METHODS 

Experiments were conducted on a 25 m x 300 m field of Volga loam 
(fine-loamy over sandy or sandy-skeletal, mixed (calcareous), frigid, 
Cumullc Haplaquoll) at the South Dakota State University Agricultural 
Engineering Research Farm located 8 km south of Brookings. South Dakota. 
Larker barley ( Hordeum vulaare L. ) was planted In the field at 15-cm 
row spacing (north-south rows) at a population of 2.5 million plants ha'\ 
The barley was not Irrigated. Surface roughness of the soil was minimal. 

Surface soil temperatures (approximately 1 mm below the soil surface) 
were measured with copper-constantan thermocouples at two locations (A 
and B) within the field. For each location, three thermocouples were 
wired In parallel to obtain an average measurement of shaded and sunlit 
soil which approximated surface temperature. For bare soil, thermocouple 
measurements were within 1.0 C of radiometric surface temperatures 
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(corrtcttd for enilssivlty and sky Irradlanct) mtasured with a procislon 

y 

radiation thtnnointtor (Modtl PRT-5, Bamos Englnttr^ng Company). Com- 
posits ttmporatures consisting of radiance contributions from the soil 
surface and the barley were measured at both locations at 1330 Local 
Standard Time (LST) on clear days at a vertical position (zerc degree 
look angle measured from nadir) Z m above the canopy. The temperature 
resolution of the 20^ field of view PRT was tO.5 C In the 8-14 um wave- 
length Interval. Canopy temperatures were measured with the PR' at a 
height of 1 m above the canopy and a look angle of 30^ from the horizon- 
tal (Millard et al., 1980) pointing to the east and the west (perpendicu- 
lar to row direction). At that angle and direction, radiance contribu- 
tions from the soil were minimized (Hatfield, 1979). Canopy temperatures 
were corrected for emissivlty and sky Irradlance using a canopy emissivlty 
of 0.98. 

Emissivltles of the canopy at full cover were measured using a pro- 
cedure similar to that described by Fuchs and Tanner (1966). We used a 
painted aluminum plate with an emissivlty of 0.52 rather than an anodized 
plate to determine sky Irradlance (81 ad and Rosenberg, 1975). Soil 
emissivltles were measured on a bare soil plot adjacent to the barley 
field. 

Soil water contents (0 to 4-cm layer) for each location were deter- 
mined gravimetrical ly on soil samples collected at the time of the temp- 
erature measurements. Percent cover was determined using 35 inn color 
Infrared slides of ths canopy (photographed from a vertical position 
approximately 1 m above the canopy) projected on a random dot grid. 

^ Mention of trade name does not Imply endorsement of a particular 
product or company. 
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Figurt 1 shows stasonal trtnds In ptrctnt covtr of th# barlty canopy. 

RESULTS AND DISCUSSION 

In the discussion that follows, composite temperature refers to 
apparent temperatures measured by the nadi reviewing PRT-5. Canopy 
temperature refers to temperature measured by the PRT-5 at a 30^ angle 
from the horizontal. 

During the Investigation, composite temperatures were 0.5 to 11.5 C 
higher and surface soil temperatures 1.5 to 20 C higher than canopy 
temperatures (Figure 2). As expected, differences between composite and 
canopy temperature decreased as canopy cover increased and less emitted 
radiation from the warm soil background was detected by the radiometer. 
The correlation between composite and canopy temperature w^.s nonsignifi- 
cant (r ■ 0.41). 

Millard et a1. (1980) found that errors from assuming nadir- viewing 
thermal scanner measurements represented actual canopy temperature were 
a linear function of canopy cover. We found a highly significant linear 
relationship (r‘ ■ 0.52) between the composite-canopy temperature dif- 
ference and percent cover (Figure 3). However, the considerable scatter 
In our data suggests that It may not be possible to assess errors in 
determining canopy temperature using only canopy cover information as 
Millard et a1. (1980) suggested. 

We assumed the longwave radiation flux from a canopy and the soil 
background could be approximated by the relationship 

R - * 0-fc)(l-Ss)B* [1] 

where R(W m ) Is longwave flux, f. Is percent cover of the canopy 
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expressed as a fraction, is canopy emissivity, is soil emissivity. 

T- (K) is canopy temperature. T (K) is surface soil temperature, a 
c ® 

(5.67 X 10’® W m’^K’^) is the Stefan-Boltzmann constant, and B* (W m’^) 
is longwave sky irradiance. The first two terms on the right-hand side 
of equation [1] represent longwave radiation emitted from the canopy and 
exposed soil background, respectively. The last two terms represent long- 
wave sky irradiance between the canopy and the soil is Ignored in equation 
[1]. Equation [1] also does not partition fractions of shaded and sunlit 
leaves, or fractions of exposed soil background which are shaded and sun- 
lit. Canopy temperature can be expressed by rearranging equation [1] to 
give 

Vc” 

We compared observed values of T with values predicted using equa- 
tion [2] and measured values of f^, T^ and B* (Fig. 4). R was calculated 

4 

from measurements of composite temperature using the relationship R " ®'*’comp 

where T is composite temperature. A measured value of 0.98 was used 
comp 

for e^. Soil emissivity varied with water content as shown in Figure 5. 

Linear regression analysis of predicted versus observed canopy temperature 

2 

yielded a slope of 1.04, an intercept of -0.53, and a r of 0.88. Differ- 
ences of observed from predicted values ranged from -1184 to +2.50 C. 

The prediction accuracy of equation [2] was independent of canopy cover. 

The correlation between the predicted minus observed canopy temperature 
difference and percent cover was 0.26 (non- significant). 
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OBSERVED CANOPY TEMPERATURE (*C) 


Figure 4. Comparison of predicted and observed canopy temperatures 
Canopy temperatures were predicted using equation (2). 




VOLUMETRIC WATER CONTENT (X) 


Figure 5 


Relationship between measured soil emissivlty and volumetric 
water content In the 0-4 cm layer. 
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Many Investigators have discussed the Importance of correcting 
radiometric data for emlsslvity variations. Barthollc et al. (1972) re- 
ported temperature errors ranging from 1.9 C for bare, dry soil to 0.8 
for cotton which arose from assuming an emissivlty of 1. Jackson c»t il. 
(1977) reported a nearly constant error of 1.7 C for wheat temperatui^ 
by not correcting for emissivlty. Similarly, Sutherland and Barthollc 
(1977) found that assuming an emissivlty of 1 produced errors on the 
order of 1.0 C for complete canopies. 

When emissivltles of 1 were used for the soil and canopy In equation 
[2], predicted canopy temperatures ranges from 6.^3 C lower to 1.70 C 
higher than observed values. Linear regression analysis of predicted ver- 
sus observed canopy temperatures yielded a dope of 1.14, an Intercept 
of -5.08, and a r^ of 0.76. Prediction accuracy varied with canopy cover 
as shown In Figure 6 with largest errors occurring at low percent cover 
when radiance contributions from the soil were at a maximum. The magni- 
tude of the error from assuming emissivltles of 1 depends not only on 
canopy cover, but also on soil type and water content. Soil emissivltles 
ranging from 0.90 for dry sand to 0.99 for loamy soils have been reported 
(Sellers, 1972; Sutherland and Barthollc, 1977; Taylor, 1979). 

Prediction accuracy when measured emissivltles were used In equation 

[2] but the B* terms were neglected also changed ,ith canopy cover (Fig. 6). 

.2 

The sum of the reflected B* components ranged from 13.2 W.m at 23^ cover 
_2 

to 5.6 W m at 90 % cover. Differences of observed from predicted canopy 
temperatures ranged from 0.8 to 10.7 C. Regression analysis of predicted 

versus observed canopy temperatures gave a slope of 0.66, an Intercept of 

2 

7.74 and a r of 0.66. 
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Figuie 6. Predicted minus observed canopy temperatures as a function 
of percent cover when values of 1 were used for and 
(circles); and when measured values of and e. were used, 
but sky irradlance terms were neglected ^(triangles). 
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This study has shown that canopy temporatures at Incomplete cover 
can be estimated using nad1r> viewing radiometers If appropriate con« 
sideratlons are given to soil background radiance, emissivlty and 
sky Irradlance. Techniques are available for estimating some of the In- 
puts to the canopy temperature model. Canopy cover can potentially be 
estimated from remote sensing data (Heilman et al., 1977; Tucker et al.. 
1978; Jackson et al., 1979). Emissivlty data can be obtained from the 
literature for a wide range of soil types. Sky Irradlance can usually 
be estimated from prevailing sky conditions (Soer, 1930). However, under 
certain conditions, sky Irradlance can be h1gh'*y variable and may require 
direct measurement (Conaway and van Bavel, 1967). Estimating the radiance 
contribution from the soil background remains a difficult problem. Models 
have been developed for estimating surface and near surface soil tempera- 
ture (Behroozl-Lar et al., 1975; Pratt and Elyett, 1979; Meyer et al., 

1975; van Bavei and Hlllel, 1976) and they can potentially be extended to 
crop canopies. 
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ABSTRACT 

Previous Investigations of thermal Infrared techniques using 
remote sensors (thermography) for estimating soil water content have 
been limited primarily to bare soil. Ground-based and aircraft 
Investigations were conducted to evaluate the potential for extending 
the thermography to approach to developing crop canopies. A signifi- 
cant exponential relationship was found between the volumetric soil 
water content In the 0-4 cm soil layer and the diurnal difference 
between surface soil tmnperature measured at 0230 and 1330 LST 
(satellite overpass times of NASA's Heat Capacity Mapping Mission - 
HCfVI). Surface soil temperatures were estimated using minimum air 
temperature* percent cover of the canopy and remote measurements of 
canopy temperature. Results of the Investigation demonstrated that 
thermography can potentially be used to estimate soil temperature 
and soli moisture throughout a complete growing season for a number 
of different crops and soils. 

INTRODUCTION 

Remotely sensed surface temperatures have been Investigated for 
estimating soil water content (Idso et al., 1975; Idso and Ehler, 1976; 
Schmugge et al., 1978). Soil water contents have been related to 
differences between the dally maximum and minimum soil or crop 
temperatures. The Investigations generally have been limited to bare 
soils or fully developed crop canopies because of difficulties in 



Interprttinq thtrmal data at Iass than full covtr whan significant 
amlttance contributions from both soil and vagatatlon occur. Tht 
ability to dtriva useful Information from remott temperature measure- 
ments for conditions other than bare soil or fully developed canopies 
would greatly expand the usefulness of the remote sensing techniques. 

Investigators have shown that even at full cover, thermal 
emittance from the soil surface can affect remote temperature 
measurements of crop canopies (Blad and Rosenberg, 1976). Thus, 
surface soil temperatures can potentially be estimated from remote 
measurements of land surface emittance where a crop canopy Is the 
primary source of radiation. 

We conducted a ground based and aircraft Investigation to evaluate 
the potential for estimating soil surface temperature and soil 
moisture from measurements of total area emittance at various stages 
of crop canopy development. The Investigation was conducted to 
examine data collected durliig times of the diurnal temperature cycle 
corresponding to data collection by NASA's Heat Capacity Mapping 
Mission (HCW), launched in April 1978. The satellite, which carries 
a two-channel radiometer (0. 5-1.1 and 10.5-12.5 ym) In a sun-synchronous 
orbit, collects data at midlatitudes at approximately 0230 and 1330 
1ST during the diurnal cycle with repeat coverage of 5 or 16 days 
dependings on latitude. 

MATERIALS AND METHODS 

Plot Study 

Experiments were conducted on a 25 m x 300 m field of Volga loam 
(fine, loamy over sandy, mixed (calcar <ous} , frigid, Cumullc HaplaquoH) 
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at th« South Dakota Statt Univtrsity Agricultural Enginotring 
Rtscarch Farm locattd 8 km south of Brookings. South Dakota. Larktr 
barlty ( Hordtum vuloaro L. ) was planttd In tht fltld at 15 cm row 
spacings (north-south rows) and a population of 2.5 million plants 
ha’^ Rainfall In tho Brookings aroa avtragts 558 mm ytar'\ No 
supplwnontal wattr was appiltd to tht barlty. Surfact roughntss of the 
soil was minimal. 

Surfact soil ttmptratures (1 mm btlow tht soil surfact) were 
measured with copper-constantan thermocouples at two locations 
(A and B) within the field. For each location, three thermocouples 
were wired In parallel to obtain an average measurement which 
approximated surface temperature. Apparent canopy temperatures consisting 
of emittance contributions from the soil surface and the barley (shaded 
and sunlit leaves) were measured with a portable 1r' ared radiometer 
(Model PRT-5, Barnes Engineering Co.) at a vertical position (zero 
degree look angle measured from nadir) at a height of 2 m above the 
canopy. The temperature resolution of the 20° field of view PRT-5 
was to. 5° C in the 8-14 urn wavelength interval. Apparent crop 
temperatures were measured with the PRT-5 at a height of 1 m above the 
canopy and a look angle of “SO® to minimize emittance contributions 
from the soil. Temperatures were measured at 0230 and 1330 1ST. 

The temperatures measured with the PRT-5 were not corrected for 
emissivity. Emissivities. determined using a procedure similar to that 
described by Fuchs and Tanner (1966). ranged from 0.96 for bare, dry 
soil to 0.98 for the fully developed barley canopy. For the range of 
temperatures and percent cover encountered, the maximum error from not 
correcting for emissivity was 1.5° C. 





1 . 


Seasonal variations in leaf area index (A) and percent cover 
(B) of the barley canopy. Jointing and heading occurred on 
16 June and 19 July, respectively. 
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Soil water contents (0-4 cm layer) for each location were determined 
gravimetrically on soil samples collected at the time of the temperature 
measurements. The average of soil water contents measured at 0230 and 
1330 LST was used to represent the 24 h average. Jackson et al. 

(1976) reported that the average of the dally maximum and **i1n1mum 
water content closely approximated the 2d h average. 

Temperature and soil water content measurements were Initiated 
when the canopy cover reached 30%. Data were collected for 22 dates 
during the 45-day Investigation. 

Plant samples for determining leaf area Index (leaf area/soll 
area) were taken every 5-7 days. Leaf areas (green leaves only) 
were measured with an optical planimeter (Lambda Instrument Corp.). 
Percent cover was determined using 35 mm color Infrared slides of the 
canopy (photographed from a vertical position ”1 m above the canopy) 
projected on a random dot grid. Dally values of lead area Index (LAI) 
and percent cover were estimated from graphs of observed LAI and 
percent cover versus date (Fig. 1). We did not estimate percentages 
of shaded and sunlit leaves, or percentages of shaded and sunlit soil. 
Maximum and minimum air temperatures were obtained from the Brookings 
National Weather Service Station ('Is km from the research site). 

All data were subjected to regression analyses. 

Aircraft Study 

Apparent canopy temperatures of corn, soybean, millet and pasture 
were collected along a 24 km flight line northwest of Brookings by a 
quantitative thermal scanner (Daedalus Enterprises, Inc., Ann Arbor, 
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Michigan) flown in the Remote Sensing Institute's twin-engine Beech- 
craft at an altitude of 3650 m above ground level. The 1.6 mrad 
scanner detector has a neat of 0.5^ C. Data wer uVlected at 1330 
and 0230 LST on 5 and 6 September 1978. Scanner data were not corrected 
for atmospheric attenuation of emissivlty variations. Sky conditions 
were clea* for all flights. Errors from neglecting emissivlty 
variations and atmospheric effects were <1° C. 

Soil water contents (0-4 cm layer) were gravl metrically sampled 
In each of the fields at the time of the aircraft overflights. Percent 
cover was estimated using the same procedure used In the plot study. 

Data from the aircraft study were used to test the predictive equations 
developed from the plot study on barley. 

RESULTS AND DISCUSSION 

Soil Water Content Versus Temperature Relationships 

The amplitude of the diurnal soil surface temperature wave is a 
function of thermal Inertia and meteorological factors (solar 
radiation, air temperature, humidity, etc.) Thermal Inertia - an 
Indication of a soil's resistance to temperature change - Is defined 
as pcx^^^, where p Is density, c specific heat and x thermal 
conductivity. Since p, c and x of a soil Increase as soil water 
content Increases, the resulting amplitude of the diurnal temperature 
wave decreases. 
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Fig. 2. Relationship of the different (aT ) between soil surface 
temperatures measured at 1330 and^0230 lST and the 
average 24 h volumetric soil water content (SWC) In the 0 
to 4 cm layer of the profile. 
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When the soil surface 1s wet, evaporation Is a major factor con- 
trolling surface heat loss. After the surface layer dries and the 
soil water supply cannot meet the evaporative demand, surface heat 
loss Is by conductive transfer (soil heat flux) and is largely 
Influenced by thermal Inertia. Nocturnal cooling Is highly related 
to thermal Inertia. Thus, the diurnal surface temperature range car 
be an Indication of soil water content. Idso et al. (1975) found a 
linear relationship between the diurnal range of surface soil 
temperatures and soil water content In the 0-4 cm layer of soil, and 
reported that the temperature versus water content relationship was a 
function of soil type. However, they also found that If soil water 
content was expressed In units of pressure potential, this 
dependence was minimal. 

Vegetation cover alters the solar radiation at the soil surface 
and thus affects soil evaporation and soil temperatures. Therefore, 
dynamic growth and development of vegetation would be expected to 
complicate the temperature versus water content relationship. 

Initially, we evaluated the relationship of day minus night surface 
soil temperatures (aT^) versus soil water content at various stages of 
canopy development. Leaf area Index and percent cover of the barley 
canopy ranged from 0.3 to 3.2 and 30-90%, respectively. The exponential 
equation 

s q(- 0.06SWC+3.59) 
s 

with an r^ of 0.81 and a standard deviation from regression of 2.54° C, 
was found to best represent the relationship between aT^ and the 
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average 24 h volumetric soil water content (SWC) In the 0*4 cm layer 
of the soil profile (Fig. 2). The exponential form fit the data better 
than linear (r ■ 0.71), power (4 ■ 0.76) or quadratic (r ■ 0.77) 
curves . 

Idso et al. (1976) proposed a procedure for compensating for 
environmental variability In the thermal Inertia approach by normalizing 
measurements with respect to an arbitrary standard diurnal air 
temperature variation. We found no significant Improvement In the 
4T versus SWC relationship using the same normalization procedure. 

The temperature versus water content relationship [Eq. (1)] 
applies only to Volga loam. However, Idso et al. (1975) converted 
soil water content to a pressure potential and found a mere universal 
relationship that appeared to be Independent of soil type. Schmugge 
et al. (1978) reported that In the absence of pressure potential data, 
textural differences In temperature versus water content relationships 
could be reduced by expressing soil water content as a percent of 
field capacity. The temperature versus soil water content relationships 
have limited usefulness unless soil temperatures can be estimated from 
remote measurements under all crop-cover conditions. 

Estimating Soil temperature from measurements of canopy temperature 

During the Investigation, surface-soil temperatures at 0230 LST 
were 1.1-2. 2° C higher than apparent crop temperatures (F1g. 3a). 
Differences between canopy and crop temperatures, even at full cover, 
probably were the result of significant amounts of thermal radiation 
from the soil surface being detected by the Infrared radiometer at 
0230 LST (Blad and Rosenberg, 1976). 
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At 1330 LSI, radiometric measurements of apparent canopy 
temperature were 0.5-17° C higher, and surface-soil temperatures 
1.5-20° C higher, than apparent crop temperatures (Fib. 3b). Greatest 
differences between canopy and crop temperatures at full cover 
occurred on days with high temperatures and high evaporative demand. 

On those days, some silting of leaves occurred, which exposed more of 
the soil background to Incoming solar radiation. 

Because emittance contributions from the soil surface apparently 
were detected by the Infrared radiometer, equations were developed 
from regression analyses to estimate soil temperatures from remote 
measurements of canopy temperature. For the 0230 LSI measurements, 
the resulting equation was 

^s(0230) ‘ ^^^(0230) * min 

with an of 0.78 and a standard deviation from regression of 1.31° C 
^s(0230) surface soil temperature, 'measurement 

of canopy temperature, and Is the minimum NWS air temperature 

For the 1330 LSI measurenent , the surface soil temperatures were related 
to the PRT measurements of canopy temperature and an exponential 
function of percent cover (PC). The resulting equation was 

^s(1330) ' '”*^(1330) * + 20.35 (3) 

with an r^ of 0.86 and a standard deviation from regression of 2.63° C 
where PC Is expressed as a fraction. We found no Improvement 1n 
estimating soil temperature by Including leaf area Index, solar 
radiation or maximum air temperature In the analyses. Fig. 4 compares 
predicted soil temperature with observed values. 



Fig. 3. Comparison of soil-surface (1 mm below surface), apparent 
canopy (PRT and apparent crop temperatures measured at 0230 
(A) and 1330 1ST (B). Canopy temperatures included emittance 
contributions from the crop and soil background. 
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Fig, 4. Comparison of predicted and observed values of surface 

son temperature at 0230 and 1330 LSI. Temperatures were 
predicted using Eqs. (2) and (3). 
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Measurements of canopy temperature used to derive (2) and (3) 
ranged from 13 to 22° C for (2), and from 24 to 52° C for (3). 

Percent cover ranged from 0.3 to 0.9. 

Evaluation of results 

Fig. 5 compires observed soil water contents with values predicted 
using Eqs. (1)-{3). and the aircraft thermal scanner measurements of 
apparent corn, soybean, millet and pasture canopy temperatures. Eq. 

(1) was converted to express soil water content as a percent of field 
capacity to minimize differences associated with soil texture (Schmugge 
et al., 1978). Percent canopy cover ranged from 50 to 80% for pasture 
and from 90 to 95% for corn, soybean and millet. Soil textures ranged 
from sandy loam to silty clay loam. Differences of observed from 
predicted values ranged from -24.5 to +15.3% of field capacity. 

The average differences was 1.6% of field capacity. The less accurate 
estimates of soil moisture for corn, soybean and millet were 
probably due to their high percent cover. 

CONCLUDING REMARKS 

The results of this investigation indicate that thermography for 
extimating soil water content can potentially be extended to 
developing crop canopies. The diurnal dlffere^ice between surface 
soil temperatures measured at HCMM overpass times Is correlated with 
surface soil water content. Surface soil temperatures can be 

r- 
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5. Comparison of predicted and observed values of 24 h average 
soil water content In the 0 to 4 cm layer of the soil profile. 
Predictions were made using Eqs. (1). (2) and (3), and aircraft 
thermal scanner measurements of canopy temperature. 




§st1mat0d from remote measurements of canopy temperature If minimum 
air temperature and percent cover of the canopy are known. Remote 


sensing evaluation of crop cover has been demonstrated (Heilman et 
a1.« 1977; Kanemasu et a1., 1977; Tucker. 1979) for certain species. 
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ABSTRACT 

Results of ground, aircraft, and satellite Investigations are 
presented that dentonstrate the potential for using data from NASA's 
Heat Capacity Mapping Mission (HCMM) satellite to provide Information 
on near*surface soil moisture;. The satellite, 

which carries a two-channel radiometer (0.5 to 1.1 and 10.5 to 12.5 ym) 
In a sun-synchronous orbit, collects data at approximately 0230 and 
1330 local standard time with repeat coverage of five or 16 days 
depending on latitude. Near-surface soil moisture Influences surface 
temperature through conductive heat transfer (affected by thermal 
Inertia) and evaporation. Thus, HCf^ data acquired near maximum and 
minimum periods of the diurnal temperature cycle can provide useful 
soil moisture Information. Hydrologic Interpretations of HCMM data 
are complicated by vegetation, evapotranspiratlon, topography, atmo- 
spheric absorption and other environmental variables such as solar 
radiation, temperature, wind, etc. 


KEY TERfIS: HCl^l, Thermal Inertia, Energy Balance, Soil Moisture 



INTRODUCTION 


Virtually all physical processes occurring at the earth's surface 
or in the atmosphere involve transformations or transfers of energy. 
Energy balance interactions have Important hydrologic implications since 
water and energy balances are Intimately related (evapotranspi ration 
requires a source of energy). Distribution of precipitation affects 
the thermal regime of the surface through changes in evapotranspiration 
and thermal properties of soil and vegetation. Surface temperatures 
are also Influenced by distribution and flow of shallow aquifers. 

Surface temperatures can provide information on the nature of surface 
and subsurface hydrology. However, spatial and temporal variations in 
surface temperature are difficult to evaluate on the ground. The 
spatial criterion can bo fulfilled by remote sensing from aircraft and 
satellite. Monitoring of dynamic hydrologic features, such as soil 
moisture, which requires repetitive coverage is feasible only with 
satellites. 

NASA's Heat Capacity Mapping Mission (HCMM) launched on April 26, 
1978, is the first satellite designed to evaluate remote sensor-deri ved 
temperature measurements of the earth's surface at times when the 
temperature variation is at a maximum. Thus, the HCMi, n-presents a 
potentially useful tool for hydrologic studies. 

Ground, aircraft, and satellite investigations were conducted in 
eastern South Dakota to evaluate the potential for using HCMM data to 
monitor soil moisture and depth to shallow groundwater. Many of the 
results are preliminary since investigations are still in progress. 
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Eastern South Dakota Is characterized by shallow perched water tables 
and significant spatial and temporal variations In soil moisture and 
agricultural land use. Most topographic features In the area are re- 
lated to glaciation or stream erosion. The complexities of the ground- 
water regime and land use patterns In eastern South Dakota provide a 
wide range of conditions In which HCMT4 data can be evaluated. 

HEAT CAPACITY MAPPING MISSION 

The HCMM carries a tt,''-channel radiometer (0.55 to 1.1 and 10.5 
to 12.5 pm) in a sun-synchronous orbit (orbital altitude Is 620 km). 

Spatial resolutions are 0.5 x 0.5 km at nadir for the visible channel 
and 0.6 x 0.6 km at nadir for the .hermal Infrared channel. The neat 
of the thermal channel Is O.A^K at 280“K. Swath width is 716 km. HCMM 
collects data at 2:3C a.m. and 1:30 p.m. local standard time at mid- 
latitudes with a repeat cycle of five of 16 days depending on 
latitude. 

Standard data products include visible, day IR, and night IR 
imagery (1:4,000,000 scale), and associated computer compatible tapes. 

An example of a night thermal IR Image Is shown in Fig. 1. Special 
data products Include day-night temperature difference and apparent 
thermal inertia (ATI). ATI, which has many attribute: of true thermal 
inertia, is defined at C(l-a)/aT where C is a constant related to 
latitude and solar declination, a Is apparent albedo obtained from 
daytime HCMM reflectivity measurements, and aT is the day-night 
radiometric temperature difference bserv*>d by HCMM. 



cmolNAL PAQC IS 
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An August 29, 1978, night HCMM thermal infrared image (scene 
ID A-AOl 25-08340) of portions of the upper Midwest. (Approxi 
mate scale 1:4,000,000, dark is cool). neAT = 0.4°C; IFOV = 
0.6 X 0.6 km; overpass time = 0234 local standard time. 





INTERPRETATION EXAWLES 


Diurnal variations of surface soil temperatures are principally 
related to thermal Inertia, evaporation, land use, and meteorological 
factors (solar Insolation, lir temperature, humidity, etc.). Thermal 
Inertia, an Indication cf a soil's resistance to temperature change. Is 
defined as (Cx)^ where C Is volumetric heat capacity and x Is thermal 
conductivity. Since C and x Increase as soil water content Increases, 
the resulting amplitude of the diurnal soil temperature wave decreases. 

When the soil surface Is wet, evaporation Is a major factor 
controlling surface heat loss since more energy Is partitioned into 
latent heat of vaporization and Is not available for heating the soil. 
After the surface layer dries and the soil water supply cannot tiieet the 
evaporative demand, surface temperature of a bare soil is largely 
related to thermal Inertia. Nocturnal cooling Is highly dependent 
on thermal Inertia. Thus, the amplitude of diurnal soil temperature 
variations can be an Indication of near-surface soil water content. 

Idso et al. (1975) found a linear relationship between the 
diurnal range of surface soil temperatures (bare soil) and near-surface 
soil water content, and reported that the temperature versus water 
content relationship was a function of soil type. The textural 
dependence can be minimized by expressing soil water content in units 
of pressure potenMal or as a percent of field capacity (Idso et al., 
1975; Schmugge et al . , 1978). Meteorological variability can be reduced 
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by normalizing the amplitude of the diurnal surface soil cemperature 
wave with respect to a standard diurnal a1r temperature variation 
(Idso et a1. • 1976). 

Vegetation alters solar Insolation at the soil surface and thus 
affects soil temperature. Therefore, growth and development of 
vegetation would be expected to complicate temperature versus soil 
water content relationships. Since crop canopies are the primary source 
of land surface emittance during most of the growing season In South 
Dakota, the use of HCMM data for hydrologic Investigations requires that 
vegetation be considered In the analysis. 

A ground study was conducted 1n a barley canopy planted In a x 

300 m field of Volga loam to evaluate soil temperature (measured at 
HCMM overpass times by thermocouples 1-mm below the soil surface) 
versus water content relationships at various stages of canopy develop- 
ment (Heilman and Moore, 1980). Percent cover of the developing barley 
canopy ranged from 30 to 90 percent over the 45-day study. The 

exponential equation 

, g(-0.06 SWC + 3.59) 
s 

2 

with an r of 0.81 was found to best represent the relationship 
between day minus night surface soil temperatures (aT^) and the 
average 24-hr volumetric soil water content (SWC) in the 0 to 4-cm 
layer of the soil profile (Fig. 2). 

The relationship In (1) has limited usefulness unless aT^ can be 
estimated from remote measurements under a wide range of crop-cover 
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Figure 2. Relationship of the difference (aT ) between soil surface 
temperatures measured at HCMM overpass times and the 

average 24*hr volumetric soil water content (SUC) In the 0 
to 4-cm layer of the profile. Temperatures were measured 
by thermocouple 1 mm below the surface In a field of Volga 
loam (Heilman and Moore, 1980). 
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conditions. The ground study found that apparent canopy temperatures 
measured 2 m above the canopy by a 20** FOV Infrared radiometer (ttodel 
PRT>5, Barnes Engineering Co.) at a vertical position (zero degree look 
angle measured from nadir) at HCfW overpass times contained significant 
emittance contributions from both the soil and the crop canopy through- 
out the growing season. Therefore, equations were developed from 
regression analyses of surface-soil and apparent-canopy temperatures 
to estimate surface-soil teniperature from remote measurements. For 
0230 LST measurements, the equation 

"^s(0230) “ *'’c(0230) ^a min ^ 

with an R of 0.78 was obtained where ^5(0230) surface soil 

temperature, and T^ (°C) Is the minimum air temperature obtained 

from the nearest National Weather Service station. For the 1330 LST 

measurement, surface soil temperature was related to apparent canopy 

temperature and an exponential function of percent cover (PC). The 

equation 

Ts(1330) • T^(,33o) x . 20.35 (3) 

2 

with an r of 0.86 was obtained where PC Is expressed as a fraction. 

Equations (1). (2), and (3) were tested using simulated HCMM 
data (aircraft thermal scanner data collected at an altitude of 3650 in 
A6L) collected over corn, soybean, millet, and pasture (Heilman and 
Moore, 1980). Percent canopy cover ranged from 50 to 80 percent for 
pasture and from 90 to 95 percent for’ corn, soybean, and millet. Soil 
textures ranged from sandy loam to silty clay loam. 
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Figure 3 compares observed soil water content with values predicted 
using equations (1) through (3) and simulated HCMM measurements of 
apparent canopy temperature. Equation (1) was converted to express 
soil water content as a percent field capacity to minimize textural 
differences (Schmugge et al., 1978). The average difference of observed 
from predicted values was 1.6 percent of field capacity. 

Preliminary analyses of actual HCMM data of eastern South Dakota 
Indicates that high soil moisture areas can be detected using HCrti 
thermal Imagery. In early April 1978 heavy spring runoff and Ice 
blockage caused significant flooding of alluvial areas In a portion of 
the Big Sioux River Basin In southeastern South Dakota (Fig. 4). Flood 
waters had receded by mid-May, but an area of high soil moisture (at 
or near field capacity) remained. Soil moisture In the surrounding 
upland soils was generally less. 

The high moisture area appeared cooler than surrounding areas 
on May 15 day thermal Imagery (Fig. 5). Temperature differences 
between the flood plain and surrounding areas were probably the result 
of thermal Inertia and evaporation differences associated with soil 
moisture differences. Landsat Imagery (Fig. 6) confirmed that no standing 
water was present In the fields. Adjacent illuvial areas did not appear 
different from uplands, Indicating that the anomaly was not associated 
with Inherent thermal Inertia of the soil but with a moisture difference. 
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Figure 3. Comparison of predicted and observed values of 24-hr 

average soil water content in the 0 to 4-cm layer of the 
soil profile. Predictions were made using equations (1), 
(2), and (3) and simulated HCMM measurement of canopy 
temperature (Heilman and Moore, 1980). 
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Figure 4. Landfom map of Brookings County, South Dakota, showing 
location of alluvial soils (bottomland) in the Big Sioux 
River Basin which were flooded in early April 1978 (Heilman 
and Moore, 1981 ) . 
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Figure 5. Lnlargement of May 15, 1978, day HCMM thermal infrared data 
(scene ID A-A0029-19575) showing a high soil moisture area 
(arrows) in eastern South Dakota (Heilman and Moore, 1981). 
Dark is cool . 
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Figure 6. Photographic enlargement of a May 13, 1978, Landsat color 
composite (scene ID E-21207-16083) of the same area shown 
in Figure 5. Note that no standing water is visible in 
the flood plain of the Sioux River Basin. 
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These results Indicate a potential for evaluating soil moisture 
using HCMH data. Final results from HCMM soil moisture Investigations 
currently In progress will! continue to evaluate the utility of using 
HCMM and similar data for quantifying soil moisture differences from space. 

DISCUSSION 

Although the potential for using HCNM and similar data In soil 
moisture Investigations has been demonstrated, there are limitations 
In the use of such data which must be considered. Environmental factors 
which Influence energy balance Interactions must be considered when 
using thermal data. Due to Its large heats of fusion and vaporization, 
water undergoing phase transformations acts as a heat source or sink. 
Changes In heat content will not be represented by a corresponding 
temperature change If a phase transformation occurs. Thus, conditions 
favoring high evapotranspi ration rates or dew or frost formation are 
not favorable for remote sensing of near-surface soil moisture. 

Wind patterns may obscure thermal anomalies created by soil 
moisture (Fig. 7). Topographic variations and vertical extrusions 
affect the boundary layer and thus affect sensible and latent heat 
transport. 

Atmospheric constituents (clouds, aerosols, water vapor, etc.) 
Influence surface temperature by attenuating Incoming solar radiation 
and affecting radiative cooling of the surface. Since atmospheric 
counter-radiation Is emitted by atmospheric constituents, radiative 
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Figure 7. Wind patterns on night thermal Imagery of an area south of 
and Including Sioux Falls. Wind wes from the northeast at 
a speed of 10 knots. Approximate scale 1:55,000; dark Is 
cool . 
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Figure 7. Wind pat:erns on night thermal imagery of an area south of 
and including Sioux Falls. Wind was from the northeast at 
a speed of 10 knots. Approximate scale 1:55,000; dark is 
cool . 
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cooling will be greater under a clear sky. Atmospheric components 
also affect the amount of longwave radlaticn emitted by the surface 
that Is detected by HC144 or other thermal sensors. 

Thermal remote sensing has an advantage of relating to subsurface 
properties since surface temperatures and emittances are a function of 
both surface and subsurface properties. These preliminary results 
indicate that observations at appropriate periods within the diurnal 
cycle can provide Information on soil moisture. These and otner 
preliminary results appear promising for development of interpretation 
models to advance the use of thermography. 
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ABSTRACT 

Four dates of Heat Capacity Mapping Mission (HCMM) data were analyzed 
to evaluate the utility of HCMM thermal data for estimating near><surface 
soil moisture In a complex agricultural landscape. Because of large 
spatial and temporal ground cover variations, HCMM radiometric tempera- 
tures consisted of radiance contributions from different canopies and 
their respective soil backgrounds. However, when surface soil tempera- 
tures were empirically estimated from HCMM temperatures and percent cover 
of each pixel, a highly significant correlation (r ■ 0.74**) was obtained 
between the estimated soil temperatures and near-surface soil water 
content. 
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INTRODUCTION 

Remotely sensed surface temperatures have been Investigated for 
estimating soil water content In the surface layer of soil (Idso et a1., 
1975'. Reginato et a1., 1976; Heilman and Moore, 1980). Generally the 
Investigations have been limited to bare soil or fully developed crop 
canopies because of difficulties In Interpreting data at less than full 
cover when significant radiance contributions occur from both vegetation 
and soil. Heilman and Moore (1980) found that thermal techniques could 
potentially be extern' }d to conditions of partial canopy cover. 

Investigations have also been limited to ground or aircraft studies 
since high resolution thermal data from satellites were unavailable for 
the appropriate periods of the diurnal temperature cycle. The Heat 
Capacity Mapping Mission (HCMM) satellite, launched in April 1978, 
was the fiist satellite devoted to acquiring high resolution thermal 
data at opv.imum periods of the diurnal cycle. The HCMM carried a two- 
channel radiometer (0.55 to 1.1 and 10.5 to 12.5 um) in a sjn-synchronous 
orbit at an altitude of 620 km. Spatial resolutions were 0.5 x 0.5 km 
at nadir for the visible channel and 0.6 x 0.6 km at nadir for the 
thermal channel. The NE^t of the thermal channel was 0.4 K at 280K. 

The HCW1 collected data at approximately 0230 and 1330 local standard 
time (LST) with a repeat cycle of 5 or 16 days depending on latitude. 

An investigation war. conducted to evaluate the utility of using 
HCMM dat? ‘valuate near-surface soil moisture for a complex agricul- 


tural landscape. 
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MATERIALS AND METHODS 

Th« study was conducted In an 8 x 24 km area within the Big Sioux 
River Basin In Brookings County In snuthcastern South Dakota. Surficlal 
deposits In the drainage basin are predominantly of glacial origin, and 
consist of end moraine, ground moraine, and outwash deposits (Ellis et 
al., 1969). 

The large study area was required to obtain a wide variation of soil 
textures and agricultural land use. Soils In the basin range from poorly 
drained silty clay loams In the flood plain to well drained sardy loams 
In the slightly elevated terraces. Major agricultural land use cate- 
gories are small grains, row crops, nayland, and pasture with field sizes 
generally less than 16 ha. 

Soil water contents (0 to 4 cm layer were determined gravimetrical ly 
on samples collected on HCMM overpass days in fields representative of 
soil and land use variations in the basin. Three samples were collected 
per field and averaged. Sampling occurred between 1000 and 1400 LST. 

The soil water content for each field was used to represent an entire 
HCMM pixel, although each pixel ultimately contained more than one land 
use. 

The number of sampling locations varied because of logistical prob- 
lems associated with collecting samples over a large area In a short 
period of time. Table 1 summarizes the soil moisture data collections. 

Percent cover at each location was determined using 35 mm slides of 
the canopies (photographed from a vertical position approximately 1 m 
above the canopies) projected on a random dot grid. When the canopies 
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Table 1. Nurrber of fields sampled for soil moisture on cloud-free days of HCMM 
overpass selected for analysis. Three samples were collected in each 
field. 


Date 

No. of Fields 

Land Cover 

June 5, 1978 

4 

corn, pasture 

July 13. 1978 

2 

corn 

August 8, 1978 

4 

corn, pasture 

September 4, 1978 

13 

corn, pasture, stubble 


Table 2. HCMM scenes 

analyzed In sell moisture study 


Date 

Time 

Scene I.D. 

June 5, 1978 

1330 LST 

AA0040- 19500 

July 13, 1978 

1330 LST 

AAOO 70- 19570 

August 8, 1978 

1330 LST 

AAOl 04- 19400 

September 4, 1978 

1330 LST 

AAOl 31 -19420 
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were too tall for the photographic procedu e, percent cover was estimated 
from visual Inspection. These data were used to prepare average percent 
cover curves for the growing season for each land use category. Dally 
maximum air temperatures ware obtained from the Brooking: National 
Weather Service Station. 

Four dates of daytime HC^W data representative of the growing season 
were selected for analyses (Table 2). Cloud-free day/night data were 
available only for July 13. but Insufficient soil moisture data were 
collected on that date to relate diurnal temperatun* variations to soil 
water content. 

Radlotfietrlc temperatures were extracted for each pixel containing 
a field where soil samples were collected by overlaying computer gray 
maps of HCMM d ,a with a Brookings County map containing the sampling 
sites. Radiometric tern*’ 'ratures were corrected for atmospheric e.'fects 
by comparing HCMM and ground measurements of Missouri River reservoir 
temperatures In central and southeastern South Dakota. Radiometric 
temperatures were not corrected for emissivlty variations. 

Percentage of each land use category for each pixel was determined 
using photo1nte<^pretat1on of a May 13, 1978, Landsat color composite 
(scene I.D. E-21 207-1 6803) superimposed on HCMM gray maps via a Bausch & 
Lomb Zoom Transferscope. Percentage of each land use within each pixel, 
and the average percent cover curves for each land use category were 
used to calculate a pixel percent cover for each date of HCMM data. 
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RESULTS AND DISCUSSION 

Surface soil temperature varia^^ors are a function of soil thermal 
properties* land cover* evaporation, albedo* meteorological and other 
factors (air temperature* solar radiation* wind* etc.). When the soil 
surface is v<et* evaporation is a major factor controlling heat loss. 
After the surface layer dries and the soil water ^upply cannot meet the 
evaporati**e demand* surface heat loss is largely influenced by thermal 
inertia. Thermal inertia* defined as (pcx) where p is density, c 
specific heat* and x thermal conductivity* is an indication of a soil’s 
resistance to temperature change. Since p, c, and x of a soil increase 
as soil water content increase?, the resulting amplitude of the diurnal 
surface temperature wave decreases. Thus, soil surface temperatures at 
maximum and miniwum periods of the diurnal t?mperature cycle can be an 
indication of soil water content. 

Idso et al. (1975) found line«r relationships between the diurnal 
range of surface soil temperatures (bare soil) and soil water content 
in the 0 to 4-cm layer of soil, and between the surface soil - air 
temperature difference and soil water content for the same depth inter- 
val. The temperature versus water content relationships were a function 
of soil type. However, Idso et al . also found that if soil water con- 
tent was expressed in pressure potential units* this dependence was 
minimal . 

Unfortunately* bare soil conditions seldom exist for any length of 
time in agricultural areas. Thus, remote sensing techniques must be 
developed for a wide range of ground cover conditions. HCMM rad'^ometric 
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tttfliperatures (correctad for atmosphtric effects but without consideration 
given to ground cover variations) did not correlate (r " 0.16) with soil 
water con-^^nt for the four dates analyzed. Generally, each pixel con- 
tained more than a single land use. Thus, the radiometric temperatures 
consisted of a combination of radiance contributions from different cano- 
pies and their respective soil backgrounds. Estimated percen* cover for 
the pixels ranged from 10 to 95 percent during the >tudy. 

Heilman and Moore (1980) found tnat soil surface temperatures 
beneath crop canopies were correlated with near-surface water content, 
and developed a procedure for estimating soil surface temperatures under 
crop canopies from remote measurements of a composite temperature. They 
used the euatlon 

T- » 0.79 T + 20.35 (1) 

where T^ (°C) Is predicted soil surface temperature, T. (°C) Is a compo- 
site temperature consisting of radiance contributions from the crop and 
soil background, and PC Is percent cover expressed as a fraction. 

Equation (1) was used to estimate surface soil temperatures from 
HCMM temperatures and pixel percent cover. Regression analyses were 
used to correlate T^ with soil water content (SWC). The following 
equation was 

T. - 58.35 SWC"*^-^® (2) 

with a r of 0.74 (significant at the 0.01 level) and a S of 2.1 °C 

y . X 

(Figure 1). 
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Figure 1. Relationship between predicted soil surface temperatures (T ) 
and soil water content in the 0 to 4-cm layers of the soil ^ 
profile. Predictions were made using equation (1) and hCMM 
radiometric temperatures. Soil water contents are an average 
of three measurements. 
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Soil water content was expressed as a percent of field capacity to ntlnl- 
mlze differences associated with soil texture (Schmugge et al., 1978). 

Soil surface temperatures are affected by day-to-day environmental 
variations (wind, humidity, air temperature, etc.). We attetiipted to 
compensate for some of that variability by relating differences between 
and maximum air temperature to soil water content (Idso et al., 1975). 
Table 3 lists maximum air temperatures for the four analysis dates. How- 
ever, we found no significant Improvement using that procedure. 

Results of this study Indicate that thermal data acquired from 
spacecraft can be correlated with near-surface soil moisture If consid- 
erations are given to spatial and temporal groundcover variations. Remote 
sensing of crop cover using multlspectral reflecta ice data has been 
successfully demonstrated (Heilman et al., 1977; Kanemasu et al., 1977; 
Jackson et al., 1979; Tucker et al., 1979; Hol>'‘?n et al., 1980). 

Reflective and thermal data have the adva..cage of high spatial 
resolution, but their usefulness Is lost in the presence of clouds. 
Microwave sensors have the ability to penetrate non-raining clouds. How- 
ever, the spatial resolution of passive sensors Is limited by antenna 
size, while active microwave sensing Is strongly Influenced by look 
angle and .iurface roughness (Schmugge, 1978). Thus, the combined use of 
reflective, thermal, and microwave sensors may be the logical approach 
for assessing soil moisture from satellites. Investigations of simul- 
taneous observations In the various spectral regions should be pursued 
to dCwermlne where each data set is unique and/or where the multiple 
estimates can be used to improve accura^ -s. 
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Table 3, Maximum air temperature for the four dates of analysis. 


Date 


^max (C) 


June 5 25.0 
July 1 3 26.1 
August 8 31 .7 
September 4 


31.1 
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ABSTRACT 

Results of ground, aircraft, and satellite Investigations are 
presented that demonstrate the potential for using data from NASA': 

Heat Capacity Mapping Mission (HCMM) satellite to provide information 
on perched water tables. The satellite, which carries a two-channel 
radiometer (G.5 to 1.1 and 10.5 to 12.5 urn) in a sun-synchronous 
orbit, collects data at approximately 0230 ano 1330 local standard 
time with repeat coverage of five or 16 days depending on latitude. 

Perched water tables influence surface and subsurface soil temperatures 
because of a heat sink effect created by the high heat capacity of 
water. HCMM data acquired at appropriate periods of the diurnal ano 
annual temperature cycle can provide useful information on shallow 
groundwater. Hydrologic interpretations of HCMM data are complicated 
by thermal inertia-heat sink interactions , vegetation, evapotranspiration, 
topography, atmospheric absorption and other environmental variablas 
such as solar radiation, temperature, wind, etc. 

KEY TERMS: HCMM, Thermal Inertia, Energy Balance, Soil f^isture, 


Groundwater, Perched Water Tables 
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INTRODUCTION 

This paper discusses the groundwater portion of a hydrologic 
investigation of eastern South Dakota using data from tlie Heat Capacity 
Mapping Mission (HCflM) satellite. Background material and Information 
about the HCMM are presented in the paper "Soil Moisture Applications 
of the Heat Capacity Mapping Mission" in the preceding soil moisture 
section. 


INTERPRETATION EXAMPLES 

Surface soil temperatures are controlled by meteorological factors 
and soi 1/water/ vegetation properties at depths within the diurnal 
damping deptli, and also by the ability of underlying soil material to 
store and transfer heat. For example, the high heat capacity of ground- 
water wuhin the depth of annual soil temperature variation produces 
a heat sink in summer and a heat source in winter win'ch reduces atinual 
temperature variations (Cartwright, 1968). Temperatures are influenced 
by aquifer thickness, rate of horizontal and vertical v/ater movement, 
and depth to the water table. 

Variations in groundwater depth do not s igni ficantly affect the 
amplitude of the diurnal temperature curve, but do shift tlie curve 
up or down in absolute magnitude (Huntley, 1978). Figure 1 illustrates 
the effect of depth to groundwater on subsurface soil temperatures 
measured in the Big Sioux River Basin in southeastern South i'akota. 
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A highly significant positive correlation (r = 0.68**) was found 
between 50-cm soil temperature and depths to groundwater of three 
meters or less. 

Myers and Moore (1972) and floore and Myers (1972) evaluated 
aerial thermography of tlie Sioux Basin and found that apparent thermal 
anomalies relai. 3 d to shallow groundwater could be detected during 
predawn hours in August and early September, the period of the 
maximum downward temperature gradient in South Dakota (Fig. 2). In 
addition, they found that the thickness of saturated sands and gravels 
corresponded closely to an apparent cool anomaly. During the daytime, 
they found that thermal patterns produced by differential ET rates, 
ground shadings, reflectances, and other factors masked thermal patterns 
produced by subsurface conditions (Fig. 2). 

Similar results are visible on HCMM imagery. The Big Sioux Basin 
appears cooler than surrounding areas on August night thermal imagery, 
primarily because of the lieat sink created by shallow aquifers within 
the Basin (Fig. 3). The Big Sioux Basin is not visible on day thermal 
or visible imagery because of the masking effect associated with land 
use (Figs. 4 and 5). Investigations are in progress to evaluate HC’M 
and similar data for evaluating depth to groundwater for the shallow 
water tables. 
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DISCUSSION 

An important cons i deration in the use of HCMM and similar data 


in groundwater studies is the interaction of soil moisture and 
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groundwater. An Increase 1n near-surface soil moisture produces lower 
soil temperatures during the day and higher temperatures at night, 
whereas, shallow groundwater produces lower soil temperatures throughout 
the diurnal cycle. Thus, soil moisture and groundwater affect surface 
temperature in the same direction during the day and in opposite 
directions at night during the late-summer period of maximum downward 
temperature gradients. Other limitations in the use of HCMM data are 
discussed in the previous soil moisture paper. 

Observations at appropriate periods of the diurnal and annual 
temperature cycle may reveal information on shallow water tables 
within the range of the annual damping depth (10 to 15 m in northern 
latitudes of South Dakota). These results appear promising for 
Jivelopment of interpretation models to advance the use of HCfW and 
similar data. 
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Figure 1. Relationship of 50-cm soil temperatures to depth to 
groundwater for row crops and pasture. Temperatures 
were measured In the Big Sioux River Basin during 
daylight hours on September 5-7, 1978. 
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Figure 2. Day and night thermal imagery of the Sioux Basin north of 

Sioux Falls. The night image (b) shows a broad cool pattern 
within the flood plain associated with subsurface conditions. 
Daytime (a) thermal patterns mask anomalies associated with 
subsurface conditions. The flood plain is delineated by the 
dotted line; numbers are thickness (m) of saturated sands and 
gravels. Approximate scale 1:60,000, dark is cool. (After 
Moore and Myers, 1972; Myers and Moore, 1972). 
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Figure 3. Photographic enlargement or an August 29, 1978, night HCMM 
thermal infrared image (sc?ne ID A-A0125-08340) showing the 
Big Sioux Basin, Note that the Basin appears cooler than 
surrounding areas, due primarily to the heat sink produced 
by sriallow groundwater within ‘he Basin. (Approximate scale 
1:1,000,000; dark is cool). 
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Figure 4. 


A September 4, 1978, photographic enlargement of a HCMM 
day thermal infrared image (scene ID A-AOl 31 - 1 9420) of the 
same area shown in Figure 3. Note that the Big Sioux River 
Basin is not visible because of emittance variation associated 
with land use. (Approximate scale 1:1,000,000; dark is cool). 
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Figure 5. A September 4, 1978, positive photographic enlargement of 

a HCMM day visible image of the same area shewn in Figures 3 
and 4. Note that the Big Sioux River Basin is not visible. 
(Approximate scale 1:1,000,000). 
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INTRODUCTION 

Soil temperatures are controlled not only by meteorological factors and 
thermal properties within the depth of diurnal temperature variation, but 
also by thermal properties of the underlying material. For example, Cart- 
wright (1968a, b, 1970, 1974) found that water tables within the depth of 
annual soil temperature variation create a heat sink during the summer which 
produces cooler soil temperatures throughout the diurnal temperature cycle. 

This heat sink does not affect the amplitude of the diurnal temperature 
variation (Huntley, 1978). The magnitude of the temperature anomaly associa- 
ted with shallow water tables Is dependent upon aquifer thickness, rate of 
horizontal and vertical water movement, and water table depth. 

Existence of temperature anomalies produced by shallow water tables has 
led Investigators to evaluate the potential of thermal remote sensing for 
locating shallow aquifers. Chase (1969) found that apparent cool anomalies 
on thermal Infrared Imagery corresponded with shallow groundwater. Myers 
and Moore (1972) found a correlation between predawn radiometric temperatures 
and aquifer thickness. Huntley (1978) reported that surface temperature 
anomalies related to water table depth variations could be separated from 
reflectance and thermal Inertial variations, but not from variations In 
evaporation r.^tes. 

We e'aluated the utility of using Heat Capacity Mapping Mission (HCMM) 
radiometric temperatures to estimate water table depth. The HCMM, launched 
In April 1978, carries a two-channel radiometer (0.55 - 1.1 and 10.5 - 12.5 Min) 
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In a sun-synchronou< orbit (orbital altitude Is 620 km). Spatial resolu- 
tions are 0.5 x 0.5 km for the visible channel and 0.6 x 0.6 km at nadir 
for the thermal channel. The NEaT for the thermal channel Is 0.4 K at 
280 K. Swath width Is 716 km. The HCMM collects data at approximately 
0230 and 1330 local standard time (LST) with a repeat cycle of 5 or 16 
days depending on latitude. 

MATERIALS AND METHODS 

The study was conducted In the Big Sioux River Basin In Brookings 
County In southeastern South Dakota (Figure 1). Surficlal deposits in 
the drainage basin are predominantly of glacial origin, and consist of 
end moraine, ground moraine, and outwash deposits (Ellis et al., 1969). 
Most groundwater In the basin Is obtained from shallow outwash deposits 
(within 10 m of the surface) and from sand and gravel lenses In morainal 
deposits. 

The Big Sioux River is In contact with the outwash deposits, and 
groundwater discharge forms the base flow of the river. Most of the 
aquifer recharge occurs from runoff from snowmelt and early spring rains. 
Groundwater levels in the basin usually rise from late March through May, 
and decrease from June through September. 

Soils In the basin are generally poorly drained In the flood plain 
and well drained In the slightly elevated terraces. Major agricultural 
land use categories In the basin are small grains (oats, spring wheat, 
barley), row crops (corn, soybeans), hay land and pasture. 

Water table elevations in the basin were measured In U.S. Geological 
Survey observation wells (Figure 1). Soil water contents (0 to 4-cr,i 
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Figure 1. Landform map of Brookings County, South Dakota showing the flood 
plain (bottomland) and terraces of the Big Sioux River Basin. The 
stars indicate locations of U.S. Geological Survey groundwater 
observation wells. 
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layer) at selected locations were determined gravimetrical ly on soil 
samples collected on days of HCMM overpass. 

Percent cover at several locations representative of the major land 
use categories were determined using 35 mm slides of the canopies (photo- 
graphed from a vertical position approximately 1 m above the canopies) 
projected on a random dot grid. When canopies were too tall for the 
photographic procedure, percent cover was estimated from visual Inspec- 
tion. These data were used to prepare average percent cover curves fc" 
each category. 

Radiometric temperatures from five HCMM scenes (Table 1) were ex- 
tracted for each pixel encompassing an observation well by overlaying 
computer gray maps of HCMM data with a Brookings County map :onta1n1ng 
the well locations. Radiometric temperatures were corrected for atmos- 
pheric effects by comparing HCMM and ground measurements of Missouri 
River reservoir temperatures in central and southeastern South Dakota. 
Radiometric temperatures were not corrected for emissivity variations. 


Table 1. HCMM scenes analyzed in water table study. 


Date 

Time 

Scene I.D. 

June 5, 1978 

1330 LST 

AA0040- 19500 

July 13, 1978 

1330 LST 

AA0070- 19570 

July 13. 1978 

0230 LST 

AA0078-09020 

August 8, 1978 

1330 LST 

AAO 104 -19400 

September 4, 1978 

1330 LST 

AAOl 31 -19420 
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Percentage of each land use category for each pixel containing an 
observation well was detennined using photointerpretation of a May 13, 
1978, Landsat color composite (scene I.D. E-21207-16803) superimposed on 
HCIfl computer gray maps via a Bausch & Lomb Zoom Transfer Scope. Percen- 
tage of each land use within a pixel, and ^he average percent cover curves 
for each land use category were used to calculate a percent cover for each 
pixel for every date of HCMM data analyzed. 

RESULTS AND DISCUSSION 

In an earlier paper we reported a highly significant relationship 
(r * 0.68'^*) between 50-cm soil temperatures and water table depths of 
3 m or less in the flood plain of the Big Sioux River Basin (Heilman and 
Moore, 1980). Subsequent analysis of additional temperature data indi- 
cated that the relationshif. could be extended to water tables as deep as 
5 m. Water table depths in terraces and uplands were greater than 9 m 
and did not correlate with 50-cm soil temperatures. Thus depths 
greater than 5 m were excluded in the analyses of HCMM data. 

Myers and Moore (1972) found that on daytime thermal imagery, thermal 
anomalies related to shallow water tables were overshadowed by vegetation 
differences (primarily differential evapotranspi ration rates and shading). 
Similar results were found with the HCMM data (Table 2). HCMM tempera- 
tures at 1330 LST did not correlate with water table depth, primarily 
because the temperatures measured were mainly those of vegetation, or a 
composite of vegetation and soil. 
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2 

Table 2. Coefficients of determination (r ) between HCMM radiometric 
temperatures and water table depth. 


Date 

Time 

r2 

June 5 

1330 LST 

0.02 

July 13 

1330 LST 

0.02 

July 13 

0230 LST 

0.03 

August 8 

1330 LST 

0.06 

September 4 

1330 LST 

0.02 


I^ers and Moore (1972) also reported that effects of vegetation 
were minimized at night and thus were able to obtain significant rela- 
tionships between radiometric temperature and aquifer thickness using 
predawn thermography. We did not find any significant correlation be- 
tween HCMM temperatures and water table depth for the July 13 0230 LSI 
data (Table 2), possibly because of the small variation In radiometric 
temperature (less than 2 C) within the Sioux River Basin. 

Heilman and Moore (1980) found that surface soil temperatures be- 
neath a crop canopy could be estimated from remote measurements of com- 
posite temperature using the equation 

T. » 0.79 e + 20.35 (1) 

where T (C) is surface soil temperature, T (C) Is a composite radiometric 
temperature consisting of radiance contributions from the soil and the 
crop, and PC Is present cover expressed as a fraction. Equation (1) was 
deveoped for measurements at 1330 LST. 

We used equation (1) to estimate soil surface temperatures from HCMM 
temperatures (corrected for atmospheric effects) and pixel percent cover. 
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and found linear relationships between the predicted soil temperatures 
and water table depth (Figs. 2-5). Highest correlations occurred for the 
August 8 and September 4 data. Summertime soil temperature gradients In 
South Dakota are at a maximum In August and early September. 

Although the correlations Improved In August and September, the 
slopes and Intercepts of the four relationships In Figs. 2-5 were not 
significantly different at the 0.01 level. Thus, data for the four dates 
were pooled, and the equation 

Tj « 26.90 + 1.30 0 (2) 

2 

with a r of 0.45 was obtained where D(m) Is water table depth (Fig. 6). 

Predicted soil temperatures were correlated not only with water 
table depth, but also with soil moisture (Flg, 7). Multiple regression 
analysis of the September 4 data yielded the equation 

Tj » 26.60 - 0.05 SWC + T.50 D (3) 

with a r of 0.87 where SWC(%) Is the volumetric soil water content 1n 
the 0 to 4-crr. layer. Increasing soil water content reduces the ampli” 
tude of the diurnal surface temperature variation through thermal inertia 
and evaporation effects which cannot be separated from heat sink effects 
using a single daytime measurement (Huntley, 1978). 

Results of this investigation demonstrate a potential for using 
satellite thermography to detect regions of shallow water tables and 
estimate water table depth if appropriate considerations are given to 
the effect of vegetation on the surface thermal regime. However, tech- 
niques for separating water table influences from those of soil moisture 
must be developed before satellite thermography can be a useful tool 
for groundwater studies. 
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DEPTH TO WATER TABLE (m) 


Figure 2. Relationship between predicted soil surface temperature and water 
table depth on June 5, 1978. 



DEPTH TO WATER TARIE U) 


Figure 3. Relationship between predicted soil surface temperature and water 
table depth on July 1 3, 1978. 
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Figure 4. Relationship between predicted soil surface temperature and water 
table depth on August 8, 1978. 
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Figure 5. Relationship between predicted soil surface temperature and water 
table depth on September 4, 1978. 
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DI^TH TO WATIt TAtll (m) 


Figure 6. Relationship between predicted soil surface temperature and water 
table depth when data from the four dates were pooled. 
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Figure 7. Relationship between predicted soil surface temperature and 
volumetric soil water content in the 0 to 4-cm layer. 
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USE OF HEAT CAPACITY MAPPING MISSION (HCf^M) 

DATA IN SOIL GEOGRAPHY STUDIES 

INTRODUCTION 

Producing enough food is an Increasingly serious problem In the 
world. Production of food, directly or indirectly depends partly on 
soil. Although Inventory of soil resources Is well along In some coun- 
tries It lags In others. 

The Heat Capacity Mapping Mission (HCMM) satellite Is a new tool 
that may be used In reconnaissance soil Inventories. HCMM was launched 
In April 1978, and collects data In the visible and near Infrared 
(.5-1.1 m) and thermal Infrared (10.5-12.5 ym) regions on the spectrum 
at a spatial resolution of 0.5 x 0.5 km. At mid-latitudes the satellite 
collects data at approximately 0230 and 1330 local standard time with a 
repeat coverage of five days. The two channels on the HCMM allow both 
reflectance and thermal properties to be used in soil studies. The ob- 
jective of this study was to compare land characteristics such as eleva- 
tion, soil texture, and slope aspect of known areas on HCMM Imagery. 

BRIEF LITERATURE REVIEW 

Lattman (1963) found valley side springs showed up more clearly on 
night time aerial TIR images than on conventional aerial photography. 
This was thought to be due to differences in the night time temperatures 
of the ground and the warmer ground water coming from the springs. 

Cantrell (1964) working with surface water found that as stream 
water begins cooling the thermal energy given off seems to warm the 
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vegetation and air thus forming a heat sink. The heat sink superimposes 
the drainage net on the imagery. 

Sabins (1969) used IR imagery for structural mapping in Southern 
California. He found that flooded fields in the Imperial Valley appeared 
warm on the imagery. Damp ground gave a cold signature. 

Cannon (1973) found that predawn IR Imagery provided a detailed 
drainage map and was a good way to inventory surface water distribution. 

Offield (1975) mapped structure of the front range and adjacent 
plains of Colorado on IR images derived from a scanner on a RB 57 air- 
craft. 

Schnieder et al (1979) used enhanced night time thermal imagery and 
digital data from a NOAA polar orbiting satellite to map drainage patterns 
and landforms in North and South Dakota. The Missouri and Prairie Coteaus, 
glacial moraines and partial drainage boundaries of major rivers were 
discerned. Analysis of satellite digital thermal data for western tribu- 
taries of the Missouri River showed north-facing slopes to be warmer than 
south-facing slopes by an average of 1.5°C. 

THE STUDY AREA 

HCMM imagery of South Dakota was examined for a number of dates. 
Visible and thermal IR images (day and night) were used. Figure 1 
shows a night IR mosaic from August for an area in east-central South 
Dakota. The same scene appears in both panels. The annotation in the 
lower panel indicates elevation, soil texture and slope aspect differ- 
ences. The oper right area labeled "warm-lower elev." is the east facing 
slope of the Prairie Coteau and the lower lying Minnesota-Red River Valley. 
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Figure 1. Predawn thermal infrared HCMM image (dark is cool) for an area in 
east-central South Dakota. 
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The Prairie Coteau lying west of this separates into two areas - the 
eastern cooler area of finer textured silty clay loam soils (Kranzburg 
series) and the wanner silt loam soils {Poinsett series). These two 
areas have about the same elevation (1700 ft.) but the temperature at the 
time the image was recorded was 49°F for the cooler Kranzburg series and 
52°F for the warmer Poinsett series. West of the Prairie Coteau the 
elevation is about 1280 ft. and the temperature is 49°F. This cool 
area has more dense clay loam soil having a high component of shale and 
has shale bedrock nearer the surface than the rest of the James Valley 
to the west. Temperature in the warmer part of the James Valley lowland 
(elevation about 1290') was about 56°F. In the lower left of the scene 
the warmer northeast and east facing coteau escarpment is seen. West 
of the escarpment is the higher lying and cooler Missouri coteau. 

Figure 2 shows the portion of the mosaic west of Figure 1. The 
warmer Missouri escarpment is clearly delineated as are the cooler 
Missouri coteau and the somewhat warmer James River Lowland. 

Figure 3 coverage lies west of the area shown on Figure 2 and en- 
compasses the Black Hills on the extreme west and the plains east of the 
Black Hills. The Black Hills appear cooler on the image while the north- 
facing slopes of the east-flowing streams and rivers are warmer. This 
verifies the results found by Schneider et al that the north-facing 
slopes are warmer than crest positions or south facing slopes. 

DISCUSSION 

The 3 figures show landscapes from night IR taken in August. The 
results show that elevated areas, finer textured soils, and crest or 
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Figure 3. Predawn HCMM ttiernial infrared image of area including the Black Hills 
of South Dakota (dark is cool). 
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south- facing slopes are cooler than lower lying areas, coarser textured 
soils, and north-facing slopes. Scenes examined for June 21, 1978, and May 14. 
1978, show the same relationships. 

Perhaps the most significant relationship for soil geography Is the 
warmer temperatures that exist on north-facing slopes. Logically the 
opposite situation would be expected since these slopes escape the direct 
rays of the sun. It appears that the lack of direct sunlight results In 
cooler temperatures during the day which in turn means less evapotrans- 
poratlon. Thus, more moisture Is available which holds heat better 
than the drier crest and scuth-facing slopes. 

No difference In soil mapping units In soil survey operations Is 
recognized on north-facing slopes from other aspect sites. One reason 
Is that prior to HCMM there was no means to measure this heat difference. 

Yet the heat difference must be a significant factor in. vegetative growth. 

Since most of western South Dakota is rangeland the difference in grass 
production and stocking rate probably is the principal kind of land use 
affected. 

Elevation differences of as much as 1000 feet occur between the 
northern part of the Prairie Coteau and the Red River-Minnesota River 
Lowland. Elevation differences are somewhat lower than this (about 
600-700 feet) between the northern part of the Prairie Coteau and the 
James River Lowland. Separate soil series have been mapped on the 
northern Prairie Coteau and the lowlands on either side due to soil 
parent material differences (texture and mineralogy differences). 

Although it was realized that temperature differences also occur 


among these areas there has been no method for showing this spatially 
until the Introduction of thermal Imagery. 
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Soil texture differences exist between the cooler eastern side 
of the Prairie Coteau and the wanner western side. The ea.«tern area Is 
covered by older glacial materials ard has a mature landscape where 
streams drain off excess water. These are the Kranzburg soils. The 
western area Is characterized by younger glacial deposits In an Irregu- 
lar Inmature landscape. Here precipitation tends to remain In the 
area since few streams exist and water drains to a nearby marsh or lake. 
This Is the Poinsett soil area. The HCMM Image clearly shows the high 
concentration of surface water present In the western iralrle Coteau 
area. In fact the warmer temperatures of this area of coarser textured 
soils may be due In part to the significant percentage of surface water 
holding heat through the night over the entire area. 
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ABSTRACT 

Experimental and theoretical Investigations were carried out relating to the 
detection and mapping of near-surface groundwater by the use of remote -sensed 
thermal emittance data (thermography). Soil temperature profiles, thermal 
emittance, soil moisture, depth to grounc^ater table, and other pertinent data 
were collected for test plots under fallow and crop cover conditions. 

Data were collected simultaneously for two plots having similar surface 
conditions except that one of the plots was Irrigated to create a difference In 
soil moisture profile. Calculations of surface temperature differences as a 
function of time for these soil moisture conditions were made utilizing a finite- 
difference model. The functional form of the theoretical ten^ei^ature 
difference with time was shown to be very close to the experimental apparent 
temperature difference for both crop cover and bare soil conditions. This 
result strongly suggests that a technique can be developed by which the effects 
of near-surface soil moisture can be separated from the ♦'Otal thermal emittance 
data by subtracting the data component having this fun^ onal form from total 
thermal emittance. This component could then be used to calculate soil moisture 
differences for a group of chosen slies. If a soil moisture profile Is measured 
at one site, the soil moisture profile could be calculated for other sites. 

A series of model calculations were carried out to simulate the effect of a 
near-surface water table at one site and a similar site without such a water 
table. These calculations show the effect of a water table is to give rise to 
a temperature difference at the surface which Is nearly constant through time over 
a current cycle. The presence of a water table thus may be detected from a 
component of thermal emittance difference which remains constant during the 
diurnal cycle when compared to a reference site. 

The model was tested utilizing widely spaced sites resolved by HCMM. 
Theoretical surface temperature differences calculated using the model agreed 
well with apparent temperature differences measured by HCMM. 
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INTRODUCTION 

Recent launchings of satellites with thermal Infrared Imaging sensors 
together with those planned In the near future hold great promise for appli- 
cation of thermal emittance data as a tool for resource management and development. 
These earth resources satellites allow time-sequential monitoring of land surface 
emittance over large areas of the earth at relatively low cost. They allow data 
to be readily available on a routine basis for use by the resource specialist 
In making management decisions. 

The potential use of thermography for monitoring groundwater uses remote 
sensing measurements of thermal emittance to estimate surface ten^)erature. Any 
factor which causes a variation In surface temperature may thus be measured by 
thermography. Near-surface groundwater Is such a factor. Its presence causes 
large changes In the specific heat and thermal conductivity of the soil. Phase 
transformations of water during evaporation or freezing also have large thermal 
effects on the energy budget of land surface and thus affects the land surface 
temperature. 

Complications with this method arise because soil temperature and surface 
emittance depend on a multitude of physical factors. Plant growth, aspect of 
slope, water table, wind velocity and other variables alter soil temperatures 
and thermal emittance In addition to variations associated with differing 
water table depths. Th js Isolation of emittance variations caused by the 
presence of soil moisture alone is very difficult. Therefore, models 
describing emittance variations associated with various physical features must 
be developed to Isolate their effects and to understand their Interdependence, 

This may allow one to compensate for their effect during data analysis or to 
schedule the collection of data when their effect 1s small. 

Another complication arises when the thermal infrared (TIR) Image Is obtained 
from satellite-borne sensors. The Image includes components of radiation 
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emitted and reflected by the surface modified by absorption and emission from 
the Intervening atmosphere. To determine exact values of surface temperature, 
corrections In the measured signal must be made which depend very heavily on 
atmospheric conditions and thus change with time. Corrections must also be made 
for surface emissivlty and reflection which also change with surface conditions. 
Thus considerable difficulty Is Involved In converting a satellite Image Into 
a surface temperature map. Temperature differences between two points on the 
earth nay be much easier to obtain with reasonable accuracy from TIR satellite 
data than to obtain exact temperatures. For example, absorption by the 
atmosphere will decrease the apparent temperature of two points but the 
apparent ten^erature difference betwjien the two points will remain nearly 
constant If the absorption Is similar over both points. Also emissivlty and 
reflectivity differences can also be minimized by making the comparisons between 
points which have the same plant cover such as two wheat fields. 

The emphasis of the effort was on relating these surface temperature 
differences to the variation In groundwater presence. Particular attention was 
placed on the use of the variations in temperature differences during the 
diurnal cycle to separate effects of groundwater tables in the top 50 cm of 
soil from that at larger depths. A technique of this type would be particularly 
applicable to satellites such as HCMM which allow more than one apparent 
temperature measurement during a single diurnal cycle. 

The technique envisioned for mapping groun.’/ater ove»" large areas using data 
collected by a sacelllte such as HCMM would be as follows: 

1. A reference site would be selected where groundwater would be inonitored 
(water table and near-surface soil moisture) on a continuing basis. 

2. The proposed model would be used to calculate ^he difference in near- 
surface soil moisture and water table depths between this site and a 
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second from the difference in apparent temperatures acquired by remote 
sensors during a diurnal cycle. 

3. This procedure can then be repeated using any group of sites and thus 
near-surface soil moisture and water table maps may be constructed. 

The general objective of this phase of the project was to develop and test 
a model which could je used as described In step 2 of the above procedure. 
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BACKGROUND LITERATURE 

The stud> was to develop a model useful for Interpreting remote sensing data 
for mapping of near-surface soil moisture and measurement of water table depth. 

TIR images have been used to locate springs and wells (i^ers and Moore, l'.)72) 
by the use of preoawn images. This together with a former study by Myers 'and 
Heilman (1969) showed that predawn linages exhibited a higher surface temperature 
for bare soil with higher moisture content in the top 50 cm. Myers and Moore 
(1972) evaluated the use of airborne thermography for mapping shallow aquifers 
using emittance patterns of predawn thermography. They obtained statistically 
significant results for predicting the thickness of the saturated sands and 
gravels for an August (maximum annual downward temperature gradient) predawn 
flight over shallow aquifers in eastern South Dakota. In a further study. 

Moore and Myers (1972) illustrated the thermal response to climatic variables 
for diurnal and seasonal thermograp'iy . Lard use, soil moisture, and other sources 
of thermal differences were easily observed for daytime thermography with their 
effects diminishing for predawn thermography . They concluded predawn August 
data were the most useful for identifying shallow aquifers in South Dakota. 

Several investigators h^ve studied the relationship between thermal emittance 
measured from aircraft altitudes and soil temperatures. Schmugge (1978) and 
Reginato (1976) have shown agreement between such TIR temperatures and those 
measured by thermocouples in contact with the soils. A study by Tunheim (1977) 
found a positive co**relation between aircraft TIR imagery and soil temperature 
fields caused by near-surface water tables associated with saline seeps. Results 
of this project showed the need for modifying the existing model to include 
the effects of near-surface soil moisture. 

The first evaluation of satellite thermography as an Indicator of soil 
moisture was performed by Moore, et al (1975). Analysis of SKYLAB data showed 
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a positive correlation between soil moisture and thermal emittance. It was 
concluded that thermal data from satellite altitudes had good potential for use 
in monitoring soil moisture and for irrigation scheduling. 

Quantitative estimates of soil moisture using thermal emittance data, 
however, require a model relating the effects of subsurface soil moisture on 
the surface temperature. Nu such model has yet been developed, although 
several similar types of models exist. A model proposed by Kahle, et al (1975) 
relates the change in land surface temperature during the diurnal cycle to the 
thermal inertia of subsurface geological materials. This model, however, does 
not allow for effects of groundwater, soil moisture , evapoLranspi ration, or crop 
cover. 

Two other models have been developed by Meyer (1972) for relating surface 
thermal emittance to the presence of shallow aquifers. These models use the 
assumption that a shallow aquifer would cause the soil temperature at a 50-cm depth 
to vary 1®C to 3®C from that of a non-aquifer region. The ability of this 
subsurface thermal anomaly to produce a corresponding surface thermal anomaly was 
investigated by use of these models. The first model simulated the development 
of a surface thermal anomaly during a single night and the second simulated the 
behavior of the thermal anomaly during several successive days. 

Each model considered heat transfer in two identical soil layers of 50-cm 
thickness. Since daily variations in the temperature are smdll at 50 cm 
(Cartwright, 1968; Carson, 1961) the lower boundary temperature of each soil 
layer was held constant. The subsurface thermal anomaly was presented by 
letting these fixed temperatures differ by an amount AT. 

The first model assumed a constant heat flux due to radiation. Using a 
finite integral transform, the heat transfer problem was analytically solved. 
Results predicted that a surface temperature difference ranging from 20% to 40% 
of that assumed at a 50-cm depth would develop in 9 hours. The rate of development 



ISO 


depended only on the thermal diffusivlty of the 50-cm soil profile. Values 
of diffusivlty for the calculations were chosen according to Sutton (1953). 

A result of this calculation was the prediction that the development of a 
surface thermal anomaly does not depend on the magnitude of the. heat flux 
radiated from the surface. 

The second model assumed a surface heat flux approximated by a rectified 
sine wave and a terrestrial radiation term as suggested by Smith (1966). No 
analytic solution was possible In this case and thus a finite-difference 
technique was used in a numerical solution by computer. Calculated temperature 
profiles showed good qualitative agreement with data taken by Carson (1961). 

One significant result was that a maximum value for the thermal anomaly 
would occur at 0700 hours. This result has recently received support 
experimentally for the case of ground water associated with saline seeps 
Aaron, et al (1976). 

The finite-difference model by Meyer is the one which was modified and 
applied to this project. 
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or, p*"iC}'< 

THEORETICAL MODEL 

The finite difference heat flow model developed by Meyer (1972) uses 
homogeneous soil profiles, each 50 cm In thickness. The 50-cm depth was chosen 
since dally variations in soil temperature are small at this depth (Cartwright, 
1968; Carson, 1961). The 50-cm soil profile Is divided Into 50 one-cm layers with 
50 equally spaced reference nodes as shown In Figure 1. The m reference points are 
usually referred to a- nodal points. Notice that nodal point 1 coincides with the 
upper surface of the slab at x » 0. The point m coincides with the otner boundary 
at X = 50 cm. The heat flux Into the surface x * 0 has been denoted while the 
heat flux out of the lower surface at x • L Is denoted as qj_. 




Figure 1. Assignment of Nodal Points and Heat Flux Terms for the Finite Difference 
Model . 
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A rectified sine wave is used to approximate heat flux (Smith, 1969). 

Since it is difficult to measure, it was treated as a parameter con^osed of a 
sinusoidal solar term and a blackbody radiation term. Its functional form was 
taken as: 

q^ ■ M sin ~ R (1) 

where M, the amplitude of the solar term, is the maximum solar radiation during 
the daylight hours. The variable t is the time of day measured from sunrise and 
L is the number of hours of daylight. The second term R is the terrestrial 
radiation term as suggested by Smith (1969). The radiation term results in a 
negative surface heat flux during the night as has been observed experimentally 

(Lettau and Davidson, 1957). The terrestrial radiation equation used was, 

R . o[t'' - (T')*]. (2) 

where R is the net outgoing radiation, o is the Stefan-Bol tzmann constant, T 
is the surface temperature in degrees Kelvin, and T' is the effective atmospheric 
temperature in degrees Kelvin (FI eagle, 1950). 

Consider a volume of material surrounding node n (n = 2, 3,“*, m-1) as 
shown in Figure 2. The volume of the material rurrounding node n is Aax where A 
is a unit surface area and Ax is the distance between nodal points. The amount 
of heat transferred from node n-1 to node n is denoted by q^_-| ^ and the amount 

of heat transferred from node n to node n+1 is denoted by q^ ^ + 1 . The heat 

stored within the vclume is given by 

■j n 





Figure Z. Energy Balance for Node n. 



153 


ror one dimensional heat transfer, the law of conservation of energy 
applied to node n results in the equation 

‘^n-l ,n “ ‘^n.n+l ^ ^ sn 

The rate at which heat is transferred between nodal points is written in finite 
difference form as 


^n-l,n 


■ - kA 


T -T 

in-1 

.Ax 


^n,n+l 


- kA 


T . , - T 
n-fl n 

Ax 


(4) 

(5) 


where T^_.j is the temperature of node n-1, is the temperature of node n, 

Tn+i» temperature of node n+1 , and k is the thermal conductivity of the 

material between the nodal points. If the conductivity of the volume element 
surrounding each nodal point is different, the conductivity between nodal points 
may be written as the average of the volume elements. Thus for equation (5) 


k - ^-1 

2 


(6) 


Equation (4) may then be written as 


q 1 ■ -( 


k 1 + k T_ - T„ 1 

n-1 n, . n n-1 

r~ > ~ji — 


(7) 


Similarly, equation (5) becomes 

^n.n+l ^2 ^ ^ • (8) 

The energy storage term expresses the rate at which the temperature of the 
volume changes. This term may be written in finite difference form as 

T ’ - T 
n n 


E » (pc) AAx . 
sn n At 


( 9 ) 


where At is the time increment, T^ is the temperature of node n at time t and 
T^' is the temperature of node n at time t+At. 

Substituting equations (7), (8), and (9) into equation (3) and rearranging 


terms yields 
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n n 
Ac 


2(Pc)j^(Ax)' 


I (k 


n-l * ^ * Vl> 


+ (k + k , , ) T I 
n n+1^ n+1 ' ’ 

Solving for the temperature at time t+At results in the equation 

'"n-l * -''n 


<k , + k )At 
n-i n 


n 2 (pc)„(Ax) 2 ^n-1 


+ —2 Q±i_ T 

2(pc) (Ax)^ 
n 


2(pc) (Ax)' 
n 


Now consider the transfer of heat at the surface x » 0. Figure 3 shows the 
volume element for node 1. 


( 10 ) 


( 11 ) 



q. ^ 


Ax 

ST 

r 


Figure 3. Energy Balance for Node 1 


The energy balance can be written as 


q “ q, - -f E . . 

^1,2 si 


The rate of heat transfer from node 1 to node 2 is 


■" 9 ~ "^1 


k. + 


( 12 ) 


( 13 ) 
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Since node 1 is at the surface, the volume of material surrounding node 1 is 

A V 

Y" A. The energy storage term is then 


T ' - T 

e' - (PC) A -1 


(14) 


Substituting equations (12) apd (13) into equation (14) and rearranging yields 


T • _ T 

-i -i . 


At 


A (pc) Ax 

X 


+ '‘l ^ '‘2 

(pc)j^CAx)' 


T, Tjl. 


(15) 


Solving for the new temperatu-e ' gives 


T,' 


2q At 


A(pc)^Ax 


(k +k )At (!< fk.,)At 

* ' ’'l (pc)jTi>TV 


(16) 


Finally, consider tl.v-- node at the lower boundary x * L. Figure 4 shows the volume 
element for node m. 


A 


i 

i 

h - 
1 ^sm 


^m- 1 , n 


_ i - 


N' in 






Ax 

,.x 


Figure 4. Energy Balance for Node ra. 

The energy balance equation for node m is 


q_ 1 “ Q, + E 

ni-l,n SIT. 


The rate of heat transfer from node m-1 to node m is 

k , + k T - T , 

q * _ ( ) A ^ . 

m-i,ni l Ax 


(17) 


(■ 


9 


(18) 
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A w 

Again since the volume element surrounding node m Is only A, the energy 
storage term can be written as 


E 


sm 


(pc) 




(19) 


Substituting equations (17) and (18) Into equation (19) and rearranging results 


In the equation 


m 


- T 


ra 


Ac 


k , + k 
tn-l m 

(pc) (Ax)^ m-1 

m 


- T - 


2 L 

9 

A (pc) -X 
m 


( 20 ) 


Solving for the new temperature gives 


T ' 
m 


(k 1 + k )At T , 
m~i m m-1 

(pc) (Ax)^ 
tn 


+ fl - 


(k_ 1 + k )At: 
-JilL-—-..™ ! T 

(pc) '^x)T~ •* ^ 

m 


(pc) .\Ax 
in 


( 21 ) 


The finite difference equations have now been derived. These are equations 
(11), (16), and (21). To solve a heat transfer problem, the Initial temperature 
of each of the m nodal points must be specified. This is identical to the 
specification of an Initial condition for an analytically solved problem. To 
calculate the new temperature at time At, the heat flux terms and q^^ must be 

specified. Equations (16) and (21) can be used to determine the new boundary 
temperatures. The new temperature of each of the Interior nodal points can be 
determined by solving equations (11) for each node. The resuHant temperatures 
obtained for the m nodal points can be used to cilculate the temperature at time 
2At. The Interaticn process Is continued to obtain the temperature at any desired 
future f*me. 

Choice of values for ax and At deperids on the thermal properties of the 
soil considered and the thickness of the scil layer. For the 50-cm layer 
considered and the thennal properties of soil used, the values At = 60 seconds 
and AX = 1 cm were found to be sufficient. 
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The model was modified during this study by considering two soil profiles 
having different soil moistures but which were Identical In other respects. When 
the percent soil moisture (ev^) and percent soil solids (&v ) are known, the 
percent air (aeration porosity Ea) can be found. With these values, the heat 
capacity and conductivity can be calculated In the following manner. 

C « 0v C + 0V..C. + E,C , (22) 

w w s s a a’ 

where C^, and C. are the heat capacities of water, soil and air, 
respecti vely. 

The values used for heat capacities are: 

- 1.00 cal/cm^/*C 
Cg - 0.48 cal/cm^/“C 
C, • 0.00030 cal/cm^/®C 

d 

Since C. Is a small part of the heat capacity. It Is neglected In model 

d 

calculations. The original model was only applicable to a homogeneous soil layer 
of 50-cm, thus It was modified so each 1-cm soil layer could have a different 
moisture value. If the moisture of the soil profile varies with depth, variations 
.n heat capacities and thermal conductivities occur. To adapt to these non- 
homogeneous conditions the model was modified to accept experimental soil moisture 
values at depths of 1, 8, 25, and 42 cm. Values are then calculated by the model 
by interpolation and extrapolation over the rest of the 50-cm profile. Thermal 
conductivity for each soil volume Is calculated by the method developed by 
DeVries (1963). This method generates an apparent thermal conductivity which 
approximates heat transfer due to mass movement of water, phase change? of water, 
convection, and conduction. 

This equation Is given by: 

^■^1 Xj Xj , 

X^ 


A 


(23) 


]58 


i 

r 


where x Is the apparent thermal conductivity of a granular material; x^ Is the 
thermal conductivity of the soil's Individual components; Is the volume 
fraction occupied by each soil fraction; and Is the ratio of the average 
temperature gradient In the granules across the medium. The value of <can be 
calculated from the following question: 


• 1/3:, [1 * l)g,r’. 

0 


(24) 


The g, value Is found using an unsaturated soil using water as a continuous medium; 

- 0.333 - ^ (0.333-0.035) (25) 

where z Is the soil porosity. 

The conductivities of the various soil constituents, x, are given these 
val ues : 


Xj ■ conductivity of soil ■ 0.00525 cal /cm sec “C 

x^ » conductivity of water » 0.00142 cal /cm sec ®C 

x^ » conductivity of air » 0.0000615 + 0.00196 x, cal/cm sec “C 

a w 

The finite-difference model (Figure 5) has the following Inputs: (1) soil 

heat flux, (2) soil moisture profile, (3) dry soil conductivity, (4) physical 
properties of the soil which Include (a) bulk density, p, (b) amount of soil by 
volume, (5) initial temperature profile and (6) effective air temperature. 

The effect of water table Is entered Into the model as a difference In soil 
temperature at the 50-cm depth. For example, a water table present at one site 
causes Its soil temperature at a 50-cm depth to be cooler by a constant amount 
In summer than an identical site with no near-surface water. 

Outputs from the model calculations are soil temperature profiles for the two 
sites and a surface temperature difference as a function of tlm*. 
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Figure 5. Schematic representation of the finite-difference model 
present format. 
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DA'A COLLECTION 

Data used In this phase of the project were of two types. The first type 
consisted of Intensive data collected on small soil plots over several diurnal 
cycles. The second was a series of data collections on a group of sites which 
were separated so as to allow resolution by HCMM sensors. Only the first type 
df data collection Is described here. 


Sites chosen for this study were located at the South Dakota State University 
Agricultural Engineering Farm, which Is near Brookings, South Dakota. Since soils 
vary considerably In this area, soil texture by hydrometer method, bulk density 
and porosity were analyzed throughout the profile depth of 50 cm. Results are 
shown In Table 1 and In Figures 6, 7, and 8. The percent of voluitw occupied by 
soil particles for a dry soil condition Is shown in Figure 6. Since the percent 
soil particles Increases with depth, the porosity decreases with depth. Figure 7 
shows the variation In soil components with depth. The bulk density Increases with 
depth as shown in Figure 8. 


Table 1. Physical Properties of Soil Used, 


Sample Depth 

Particle Size 
Sand Silt Clay 

Porosity (Eff) 

Bulk 
♦ i i;y 

(cm) 

(%) 

{%) 

(X) 


(g/cm^) 

0.0 - 7.6 

27.7 

61.5 

10.9 

.49 

1.36 

15.2 - 22.9 

25.0 

65.5 

9.5 

.47 

1.40 

30.5 - 38.1 

18.6 

74.3 

7.1 

.41 

1.47 

45.7 - 53.3 

15.6 

78.4 

6.0 

.39 

1.61 


2 

Each data collection site was divided into two plots, each approximately 10m . 
To prevent water movement from one plot to the other a trench was excavated to a 
depth of 100 cm and a plastic barrier buried. This barrier allowed one plot to 
remain dry while the other was irrigated to whatever soil moisture desired. 
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Figure 6. Percent soil solid by volume of soil profile as a function of depth. 
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Figure 7. Percent of the soil particles composed of sand, silt and clay as a 
function of depth. 
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Figure 8. Bulk density of soil profile as a function of depth. 
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To ineasure soil temperatures, thermocouples were Implanted In each plot at 
depths of 1 , 5, 10, 25, 50, and 100 cm. Thermal emittance (apparent surface 
temperature) was measured utilizing a Barnes PRT-5 mounted on an appartus as shown 
In Figure 9 which scanned each plot every IF minutes during data collection. Soil 
and air temperatures together with relative nw dity were collected every hour 
while solar radiation and net radiation data were collected every fifteen minutes. 

Soil moisture data were acquired by the gravimetric method with collection of 
sol.' samples an hour before solar noon (Jackson et al., 1976) to best represent the 
average moisture content. The gravimetric method of soil moisture gives a value of 


soil moisture by weight, e^, 

9 3 mass water 
®m mass dry soil 

*1 soil moisture b 
density to give the model Input or. 


(26) 


In the model soil moisture by volume Is required. Thus Is multiplied by bulk 


‘ 27 ) 

where 9 Is volumetric soil moisture. 

V 

Data were collected for several diurnal cycles for both barley and rye crop 
covers together with bare soil conditions after these canopies were removed. 

Figure 10 shows the barley crop canopy for which detailed data are used In this 
report. Data for the rye crop were qualitatively similar and will not be shown 
In detail. 
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Figure 9. 


Photograph showing the scanning apparatus used to move the 
Barnes PRT-5 across the experimental plots. 


Figure 10. Barley crop canopy present during the data collection 
on August 5 and 6, 1978. 
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RESULTS 

Data were collected for bare surface conditions for a 52-hour period starting 
at 1000 hours, August 7, 1970. Two adjacent plots were prepared as previously 
described so that one plot would have a higher soil moisture profile than the other. 
Gravimetric soil moisture measurements for these plots were made as a function of 
depth. A smoothing of the curve of field measurements was conducted using a cubic 
spline as described by Kimball (1976). 

Apparent surface temperatures of these plots are shown In Figures 11 and 12 
for the 52-hour period. The measurements were made utilizing the Barnes-PRT 5 
mounted on the scanning apparatus shown In Figure 9. The points shown are the field 
values while the continuous curve resulted from smoothing data by use of a cubic 
spline. Note the amplitude of the temperature variation during the diurnal cycle Is 
less for the higher moisture plot. This is consistent with the results reported 
previously by Idso et al (1975). 

The apparent surface temperature difference between the two plots is shown as 
a function of time In Figure 13. Values for this plot are temperature differences 
calculated from the values of the cubic spline curves of Figures 11 and 12. A 
surprising feature Is Its close similarity In functional form to the individual 
apparent temperature curves from which it was derived. This similarity Is 
particularly significant since the diurnal amplitude of the surface temperature 
has previously been related to near-surface soil moisture (Idso et. al . 1975; 

Idso and Ehler, 1976; Schmugge et. al . , 1978). Since the temperature difference 
curve shows the same functional form as curves of surface temperatures, the 
amplitude of the temperature difference should also be related to soil moisture 
differences. This type of technique for remotely measuring soil moistures would 
have the advantage of bypassing the calibration problems inherent in thermal 
emittance measurements. 
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Figure 11. Apparent temperature for bare soil surface of the dryland plot soil 
profile beginning 1000 hours, August 7, 1978 and continuing for 52 
hours thereafter. 
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Figure 12. Apparent temperature for bare soil surface of the irrigated soil profile 
beginning 1000 hours, August 7, 1978 and continuing for 52 hours 
thereafter. 
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Figure 13. Apparent temperature difference between the dry and irrigated plots 
for bare soil. Data is shown for a 52-hour period starting at 
1000 hours, August 7, 1978. 
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Results of theoretical model calculations for August 8, 1978, are shown In 
Figure 14. Inputs to the calculation are the bulk densities of the soil, the 
measured net radiation, and percent soil moisture by volume (See Appendix A). 
Comparing Figure 14 with the experimental plot of Figure 13 shows the functional 
dependence of the theoretical curve to agree very well with the experimental 
curve, particularly during the daylight hours. The magnitude of the calculated 
daylight temperature difference Is smaller with a maximum calculated temperature 
difference of about 7® C compared to a measured difference of about 10® C. A 
possible explanation of this difference Is the additional cooling of the 
Irrigated plot due to water evaporation from the surface. Increasing the amplitude 
for net radiation allows the model to simulate the daytime surface temperature 
differences very accurately. However, the calculated soil temperature profiles 
become much wanner than those measured. 

Experimental soil temperatures are compared In Figures 15-18 with theoretical 
temperature profiles calculated by the model using measured Inputs. The functional 
form of the calculated temperature profiles are very similar to the measured 
profiles. The theoretical values, however, tend to be warmer during the day and 
cooler during the night. This result also implies that evaporation from the surface 
cannot be Ignored in model calculations and must be accounted for with a parameter 
which effectively reduces the net radiation term to obtain the soil heat flux. 

Data were collected for the barley canopy shown In Figure 10 for these same 
two plots prior to the bare soil data previously discussed. These plots were 
prepared in the same general manner as the bare plots and the same types of data 
were collected. Figures 19 and 20 show the apparent surface temperatures of the 
two plots for a 32-hour period beginning at 1100 hours, August 5, 1978. Figure 21 
shows the apparent surface temperature difference obtained by subcontracting cor- 
responding apparent temperatures from the cubic spline graphs of Figures 19 and 20. 
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Figure 14. Theoretical model calculation of surface tanperature difference between the dry and irrigated 
plots for a 52-hour period starting at 1000 hours, August 7, 1978. 
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Figure 15. Soil ten.perature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (6) for August 7, 1973 at 
0400 hours. Experimental temperatures art represented by (++). 
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Figure 16. Soil temperature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
0800 hours. Experimental temperatures are represented by (++). 
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Figure 17. Soil temperature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
1600 hours. Experimental temperatures are represented by (++). 



Figure 18. Soil temperature profiles calculated by the theoretical model for 
the dry plot (A) and the irrigated plot (B) for August 7, 1978 at 
2200 hours. Experimental temperatures are represented by (++). 
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Figure 19. Apparent temperature measured over the barley canopy of the dryland 
plot. Data collection begins at 1100 hours, August 5, 1978 and 
continues for 32 hours thereafter. 
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Figure 20. Apparent temperature measured over the barley canopy of the 

irrigated plot. Data collection begins at 1100 hours, August 5, 
1978, and continues for 32 hours thereafter. 
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Figure 21. Apparent temperature difference between the dry and the irrigated 
plots with a crop canopy. Data is shown for a 32-hour period 
starting at 1100 hours August 5, 1978. 
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Comparison of Figures 19 and 20 with Figures 11 and 12 shows a drastic effect 
of the plant canopy on apparent surface temperatures. Both the functional 
dependence and the actual apparent temperature values are quite different for the 
barley canopy. However, the functional dependence for the apparent temperature 
differences as shown In Figures 13 and 21 are considerably more alike In functional 
form than the apparent surface tenperatures . During the daylight hours the barley 
plots exhibit approximately one half the temperature difference of the bare plots 
but the only difference In functional form Is a slightly slower rate of decrease 
In temperature difference late In the afternoon. The crop canopy, however, 
essentially eliminates the observed temperature difference for the nighttime hours. 
These results suggest that apparent temperature differences during the middle of 
the day may be the most likely Indicator of soil moisture differences In the case 
of a thick plant canopy. 

A theoretical calculation of the surface temperature difference for the barley 
plots Is shown In Figure 22 tor August 5, 1978. Again the calculated temperature 
difference 1s smaller than the measured apparent temperature difference; but the 
ratio of the two Is approximately the same as for the wet/dry bare soil discusaod 
previously. Since no apparent temperature difference Is observed for the barley 
canopy during the night, the model obviously Is not valid in Its present form for 
that time period. 

Several calculations were carried out to determine the dependence of the 
surface temperature difference on surface soil heat flux. In these calculations 
two plots were considered with soil properties identical to the experimental plots 
used for this study. Soil moisture by volume was assumed to be 10% In one plot and 
20% In the second. Soil heat flux values were chosen to span a range which would 
Include most experimental situations for a clear day. The maximum temperature 
difference predicted during the day was plotted as a function of the maximum soil 
heat flux. Results shown in Figure 23 display the resulting relationship. If 
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Figure 22. Theoretical model calculation of surface ten^erature difference between the dry and *>r jested plots 
for a 32-hour period starting at 1100 hours August 5, 1978. 
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Figure 23. Calculated surface temperature difference between two plots as a 
function of soil heat flux amplitude. One plot is considered to 
have a constant soil moisture profile of 10% by volume while the 
other has a constant profile of 20%. 
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further theoretical and experimental results show this relationship to be valid, 
differences In dally solar radiation which exist during the satellite overpass 
could be easily accounted for during analysis of data. 

Calculations were also earned out to determine the relationship which would 
be expected to exist between moisture difference and maximum temperature difference 
observed during the day. Values for soil heat flux and physical properties of the 
soil were again chosen to correspond to the experimental plots of this study. The 
refei-ence plot was chosen to have a soil moisture of 10% by volume and the soil 
moisture of the other was varied to a maximum of 21%. Results shown In Figure 24 
display an approximately parabolic relationship. If this proves true In further 
theoretical and experimental investigations, the development of a paractical 
technique for utilizing apparent surface temperature differences to measure soil 
moisture will be greatly simplified. 

A series of calculations were carried out to determine the feasibility of 
using dirunal surface temperature fluctuations (thermal inertia) as a measure of 
presence of water tables below the 50-cm profile. Two identical soil profiles were 
considered with the temperature at a 50-cm depth for one of these profiles cooler 
by an amount aT due to the presence of groundwater. Both profiles in these 
calculations were assumed to have the same soil heat flux at the si-face. 

The amplitudes of the surface temperature variations for a diurnal cycle were 
the same for both profiles. The surface -emperat the profile with the cooler 

50-cm temperature, however, remained cooler throughout the diurnal cycle by a 
constant amount. This constant surface temperature difference varied between 20% 
and 50% of the 59-cm difference of aT depending on soil parameters. 

These calculations suggest the variation in surface temperature differences 
during a diurnal cycle does not depend on the temperature at a 50-cm depth, but 
rather on soil parameters within the 50-cm profile and the soil heat flux. Thus 
thermal inertia would seem an unlikely candidate for determining presence of 
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Figure 24. Calculated surface temperature difference between two soil plots as a function of moisture 
difference. One plot is considered to have a fixed soil moisture profile of 10% by volume 
while the other is varied from that value-. 
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subsurface water tables unless the presence of th(*se water tables affect the soil 
moisture In the 50-cm profile. However, the presence of such a subsurface water 
table win give rise to a constant temperature difference at the surface. In order 
to separate this factor from other factors which would give rise to time dependent 
teirperature differences one will need a minimum of two observations of apparent 
surface tenperature differences during a diurnal cycle which a satellite such as 
HCMM makes available. Also data from several successive days would be helpful In 
removing other factors from the data. 

Calculations were made for the series of sites where ground truth was 
collected at the time of HCMM overpass. For these calculations average soil 
parameters for the area were used as model Inputs together with experimentally 
measured soil moistures and 50-cm soil temperatures. Soil heat flux was adjusted 
In the model so that soil temperatures calculated agreed approximately with 
experimental values. This value was then used for all calculations for the 
respective site. To account for differences In crop cover, the value of soil heat 
flux was adjusted for one of the two sites. 

Results of the calculations are shown In Table II. Calculated surface 
temperature differences agree well with the apparent surface temperature differences 
from HCMM. To use the model as a method of estimating near-surface soil moisture 
and depth to water tables, the differences In surface soil heat flux was used as 
a parameter in these calculations. 

Therefore, differences in soil heat flux must be measured, predicted by a model, 
or empirically estimated using readily available data. Since percent cover, 
differences, temperature difference at 50-cm, and apparent surface temperature 
difference data were available for this investigation, a step-wise multiple 
regression (28) was conducted to determine their use in predicting surface heat flux 
»» differences using data presented in Table II. 
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y - (7.64 + 1.56 X 2 - 1.0 X 3 - 5.2) (10^) (28) 

where y 1 s soil heat flux difference, x Is percent difference of canopy cover, 

X 2 Is percent soil moisture difference, and x^ Is 50-c«. soil temperature 
difference. The resulting equation when entering variables significant at the 
0.05 level follows with a multiple step-wise correlation coefficient as 0.835 » R. 

Therefore, up to 70% of the variation could be accounted for using measure- 
ments or model predication which could be available during model Implementation. 
Percent cover could be estimated with Landsat, apparent surface temperature 
difference between two sites estimated with HCMM or other satellites, and the 50-cm 
temperature difference estimated with a model prediction. Soil moisture and heat 
flux differences may possibly be separated by use of multiple data collections of 
surface apparent temperature differences during the diurnal cycle. The technique 
shows promise and should be explored further. 
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PROJECT CONCLUSIONS 
Conclusions of this stu(ly are: 

1. Soils with different soil moisture profiles differ In surface temperature In 
a well-defined functional manner during the diurnal cycle. This functional 
dependence Is similar to the diurnal surface temperature variations for each 
plot. 

2. A thick crop canopy destroys the water-related apparent surface temperature 
difference during the nighttime hours. 

3. The functional form of the apparent surface temperature difference measured in 
the field diurnally which was associated with a soil moisture difference Is 
changed less by the presence of a crop canopy than the functional form of the 
Individual surface temperatures. 

4. The theoretical model used in this study predicts a functional form for the 
apparent surface temperature difference very similar to that observed for the 
daylight hours. However, the magnitudes of the theoretical temperature 
differences are smaller than the experimental values for both a bare soil and 
a crop canopy. The ratio of calculated temperature difference to that 
measured is approximately the same in both cases. Multiple pass satellite 
dat? would serve a better purpose if acquired solely during daylight hours 
rather than in an orbit similar to HCMM where a night pass is included. 

5. Since the observed nighttime surface temperature differences vary considerably 
in functional form for the bare soil situation and are zero for dense crop 
canopies during the night, nighttime emittance data do not seem promising for 
use in measuring soil moisture. 

6. Model calculations predict a linear i'elationship between soil heat flux 

and the surface temperature difference arising from soil moisture variations. 

7. Model calculations predict a relationship between surface temperature 
difference and soil moisture difference. Only limited field studies were 
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conducted during this effort and further field verification should be 
pursued. 

8. Model calculations Indicate that day-minus night apparent cemperatures are 
not useful In predicting depth to the water table. However, either day or 
night temperature differences between two sites are related to differences 
in water table depths. Field data showed that where dense crop canopies 
were present, the night data differences are not observed. 

9. Model calculations of surface temperature differences agree well with 
observed apparent temperature differences from HCMM when soil heat flux Is 
used as a parameter. 

10. Initial observations are that a significant portion of the variance of soil 
heat flux can be accounted for with remote sensing observations and model 
estimates. Further Investigation Into the app»^ach is warranted. 

11. The overall results of this study reveal promise for the development of a 
method to monitor soil moisture by satellite. Using points on the curve 
comparing apparent surface temperature differences, one could calculate 
soil moisture differences for a group of chosen sites. If soil moisture 
is then measured at one site, soil moisture nay be calculated for the 


other site. 
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UPL PR0GKA.M l-OK I'lMTn-O: HEAT FLOIV SIMUI.A I ION MOHLI, 


The following program lij>ting is writ ion in a lanjuagc usc'i by 
Hewlett-Packard in the 9tJ35A mini -computoi futiiiEhcJ by the h'utcr 
Resources Institute at South Dakota State University 


C(*J Conductivity Profile A 

D(*l Conductivity Profile B 

E[*J Heat Cuj city Profile A 

F[*) Heat Capacity Profile B 

H(*] Time of Day (Hour) 

A[*] Temperature of Profile A 

B(*] Temperature of Profile B 

CfO] Term us?d by DeVries in calculation of coiHliicuivit) 

G[l] Heat capacity of water 

G[2] Heat capacity of soil 

H(l] Initial starting hour 

H(2] Ending hour 

I[0:S0] Tliermal inertia for site A 

J[ 0:501 Tliermal inertia for site B 

K[0] Conductivity of air 

K(l] Conductivity of water 

K(2j Conductivity of soil 

Mil] Initial starting minute 

M[2] Minute when calculation is to end 

Of*] Soil woisturc ff"" site A 



PM 

Q(*) 

Rf‘) 

S[‘] 

T(‘) 

U[0-S1] 

V(0:50l 

K(l] 

W[2l 

Wl3) 

v;[4] 

W[5] 

X[x] 

Y(l) 

Z( 

A 

B 

C 

n.E.F.G 

1 

J 

K 

L 

M 

N’ 

r 

Q 


0R1G<WA\- 

AS POOR OUA.LlT^ 


Soil moisture for sLto B 

Thermal diffusivity for site A 

Allocation to store old temperature for site A 

Allocation to store old temperature for site B 

Thermal diffusivity for site B 

Aeration porosity for site A 

Aeration porosity for site B 

Effective air temperature 

Amount of soil by volume at 1 cm depth 

Amount of soil by volume at 8 cm depth 

Amount of soil by volume at 24 cm dcptli 

Amount of soil by volujne at 42 cm depth 

Distance between nodal points 

Ending day 

Real data 

Soil heat flux for site A 

Soil heat flux for site B 

Time from sunrise to solar noon 

Soil bulk density at 1, 8, 24 and 42 centimeters 

Count or 

Counter 

Counter 

Day Length 

Amplitude of soil heat flux for site A 
N'uiiiber of equally spaced nodal points 
Time Let ween printouts 
Amplitude of soil heat flux for sice 3 


< 
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U IXifr./ V3i* able 

T TiPic J,..‘.crviil between calcul In re'xnuls 


X DuLiTty variable 

Z Heal iata file 


OF V’ -V- 


- fc 


3 



OMQINAI. mot 18 
Of fOOR QUALITY 


PROGRAM LISTING (IIPL) 


Operational Procedure 

Ot 70&*r0i9to 119 
*103U2 


Subroutine for Soil Hct Flux 

It "OHCNT-tMlll /60 t 21 
2t I.(l)*r2l-»r22 
3t r22-C*r23 

4t (Wm+273.\G)/l00-c24 
5: (N|0)+273.1G)/I00^r25 
6t (n(OI+273.13)/lOO*r26 
7t .000136»r24“lT27 
8t -(.000133*r25“4-r27)Tll 
9i -(.000135*c26*4-r27)Tl2 
10: if r23<«0;gto 15 
11: if r23>*L;gto 15 
12: rnd 

13: (r23* n/L) fcllTll 

14: 0*nin (c23* TT/L) +rl2*rl2 
!.5: rll ♦A; 1 12* S} re t 
'*17413 


Subroutine for Construction of Profile Plot 


16 

17 

13 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

3 : 

36 

37 

38 

39 
49 

41 

42 

43 

44 

45 

46 


"Plots": 
dcg;0«R 
sci 0, 10«3,7 
f<d 0 

C3l2 1.25,1,1,0 

pit .5,6.5,l;lbl "SOIL TRMP POR SITE A" 
pit 3.7,6.0,i j Ibl "ilOUR " 
olL 4. 7, 6. 3,1 

if HUKIO; str(H{l n*K$; Ibl "0"}colt -l,0;lbl ES,"00"}jno 2 

cnlt -l,0;otr(H(lJ )*n?;lbl r.S,"no" 

pit 5. 5, 6. 5, 1 ; Ibl "SC’L TEMP FOR SITE 13" 

pit 1,1,1 

pit 4,1,2 

nit 4,0,2 

pit 1,0,2 

pit 1,1,2 

pen 

csiz 1,1, 1,0 
foe to -50 by -5 
If I<-‘)jclt . 4, 6 H*. 1,1; Ibl I;jno 2 
olt .55,6+r".l,I;lbl I 
pit .95,5+I*.l,l;lbl "-"jnoxt I 
CSlz 1.25,1,1,99 
pit .2, 2. 6,1; Ibl "DEPTH (CM)" 
csiz 1,1, 1,0 
for I>10 to 50 by 5 


olt -.05+I*.075, ,75,l;lbl I 
olt .254-1* .075,1 , l;lbl "l";p.LKt I 
csiz 1.25,1,1,0 
pit 2,.5,1;U1 "TCMP (C)" 
if 0>R;ofs 5, 0 : 10* R ,*g to 27 
47: ofs -5,0; ret 
•31254 
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ORIGINAL PAQE IS 
OF POQR QUALITY 


Subroutine for Plot of Tlieoreticnl Profile A 
40 1 "olot\"t 

49t Pit . 25+A(iP.Ci75,S + Lm‘.l 
SO; if l»S0;pen 
Sit ret 
*213 


Subroutine for Plot, of Experimental Profile A 

S2i "reo?A": 

53i colt -.33, -.25 

54; r.lt . 25+2(J, :i [1 1/7., 1 1 * .075, 6-M U ) * . 1, 1 j Xbl 

S5t ret 

*10070 


Subroutine for Plot of Theoretical Profile B 
SG: -olotn": 

57: pit 5.25+B|Il*,a75,6*r.lll*.3 

50: if I«50rpoa 

53: rt^t • . 

*1007 « 


Subroutine for Plot of Experimental, Profile B 

60: "lealO"? 

61: colt . 33, -.25 

62: olt 5. 23*-S(J,!Ull/2,2J*.075,6-:J[Jl*.l,l;lbl 

63: ret 

•10725 


I 


C 


OWGlNAl- P*®* ® 
OF POOR 


Subroutine for Construction of Temperature Difference Plots 


64: "riff olot": 

66: dcq;0*n 
66: ec 1 OiXO^Of# 

67: csiz 1.26,1.1,770 

68: pit 4.5,6.2,l j Ibl "SURFACn TC‘»i»3RATUilC DlPrEHn JCfi' 
69: olt 9. 0,6, 1 1 Ibl "5 CM TSMHERAT UPC JirPCRCNCE" 

70: Cxd 0 

71: pit 1,1,1 

72: olt 4,1,2 

73: pit 4,G,2 

74: Dlt 1,6,2 

75: pit 1,1,2 

75 2 pen 

77: CRiz 1,1,1,270 
78: tot X«0 to 24 by 4 

79: If K10;olt .7,5. 2-I*.2i,l:lbl Ijjmo 2 

eO: nit .7,5.25-I*.21,l»lbl i 

81: pit .9,6-I*.21,l;lbl "|"»next 1 

82: C8iz 1.25,1,1,270 

83: pit .5,3.75,lilbl "HOUR" 

34: esiz 1,1,1,270 
35: for X»-5 to 10 by 2 

86. if K9}olt 2.125*1* .1875, 6.5, l»lbl X," -"jjmp 2 
37: pit 2.125+I*.1375.6.6,l;lbl I," 

88: next X 

89: cslz 1.25,1,1,0 

90: olt 2.2,6.73,1 jlbl "TEMP (C)" 

91; lino 1,2 

v2i for X 'O to 24;olt 2. 125,6-I*.21}noxt Itoen 
93: fxd 2 
94: line 

95: if R»Ci;ofa 5,0>1' *»<jgto 70 

96: o£s ■*5,0}rct 

*26454 


Subroutine for Plot of Surface Temperature Difference 

97: "olot:": 

90: line 

99: pit 2.125+U[I1*.1375,6-0(X1*.21 
100: if I»72;oen 
101: ret 
•17461 


Subroutine for Plot cf S cm Temperature Difference 

102: "olotC": 

103: lii.c 

104: Pit 7.12r>+vm*.1375,G-om*.21 

105: if l*72;nen 

106: 

• 68 “: 


ORIGSNA*- PAQE 
OF POOR QUAUrf 


Subroutine for Conductivities in Profile A 
107 » “Conu: \"j 

1031 .3T3-UU )/n-Wfxn* (.333-.0J5)-ClOl 

109 1 ,003061'i + . onl ;^*0| I )*K(01 

UOt (2/(l+(K(2)/K|l)-l)»3(01)fl/{l+(rr2)/i:m-n* (l- 2 *C(on 

111: <2/(?>(!;|O)/N|ll-l)*C(0n V(l+(Klu)/MlJ-l)Ml-2*r.lO| ) ) )/3-i*2 

112: ret 

•31043 


Subroutine for Conductivities in Profile B 

111# ft**# 

114 : . 333 -V(I)/(l-;Mxn* (. 333 -. 035 )*ClOl 
115 : . 0000515#. 00193*P( ri*K(OJ 

116 : (2/(l+{K(31/KlL)-l)*3|0nfl/{l#(K(31/Ktll-l)* (l-2*Gf01 ) ) )/3-r3 
117 : (:V(l+(KlO}/KU)-l)*Gfon+l/(l*(K(01/Kll)-l)*(l-2*f.|0) )) )/3-r2 
110: ret 
*13617 


Dimension Statements 


119: aim Ml) :501 (0:501 ,C (0:501 ,0(0:50] ,E (0:50) ,P (0:501 ,0 10 :1 05), P(0 :5C) 

120: dlM xm ,V(0:1051 ,;j(51 ,3S(41 

121: dim ‘1(0:4 ),! (0:21, K (0:3 I, C( 0:21 ,U0:501 ,Y(T; 11 ,R(0: SCI ,S( 0:51) ,f)(i.*; 105 
122: cii.m t { 5, 12 , 2) f ’ 1 , 1 ( 0; 50 ) , J( 0 : 50 ) /3( 0 : 50 1 ,7 (0 : 50) 

*10727 


) 


entering of Calculation Parameters 

123: c-nt .X!3J ,7,n.P,M(r 1 ,H (1) ,M(21 (2) , Y(1 ) ,'»v'( D 

124: ent “real data file",2;ldf X,2(*1 

*32553 


Temperature Data Jtorage 

125: l-IJ(l) j5*;i(2) }10*;:(3] ;25*:4(41 ;50«:J(5) 
126: “1*L(11 

127: £or I«2 to 50 ; L( I-l 1 -1*L( 1 1 : next I 
129: 0*J*M{3)*M(41-rl0 
129: 1!(1)*G0#:(1)^M(3) 

130: l*rl3 
*27980 


Entering of Soil Tcviiporature Profiles 

131 : ent MO) ,MIi ,A (51 ,A( 10 ),A( 25 ),A( 50 J ,b( 0 ] , 0 ( 1 ) , D( 5 ) , D (10 1 , 5 (25 ) , P. (50 ) 
*8924 
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origin;- ■ i 
OF po( ■' ■: '*•?’ 

Interpolation of Soil Tompcroturos 

132j for 1-2 to 4 > (M01-MH )/4* (l-l) 

133: (OI51-Min/4* U-lk^om -nm jnoKt I 

134: for I-S to (A(lO|-M!i I )/5MI-5)->A|5)<*NllJ 

135: (0ll0l-»I5h/5MI-U+3t5|*3m jnext I 

I36i ftr I-ll to 24>(M?.51-A(lOn/l5*|I-l)+MlOI*A(I| 

137: (QI2SI-DU0U /IS* (1-I0)i3(101 ‘nm >noxt I 

138: Cor 1-26 to 40 j (M501 -.\(251 )/25* (I -255 M (25) *M I) 

130: (3(S3)-a(25))/25' (l-2S)t3(2S)«J|l)}noxt I 

•4089 


Surface Temperature Difference 

140: A(0)-a(0)-S(Sl] 

•5957 


Determination of Mocle for Entering Heat Capacity and Thermal Conductivity 

141: J-tl*J;i£ J>I;gto 250 
*25973 


• Entering of Soil Physical Properties 

142: eat 0 (1) , D(9 ) , 0(25 ) , 0(4 2 ) , PI 1 1 , Pl» 1 , P (25) , P (4 2 ) ,W( 2) ,W(31 ,W|4 ) ,W(51 
143: ent KID ,K| 2) , K (31 , G( 1 1 , G ( 2) , D, K,F , C,M,0 ,C,t. 

•29564 


Calculatioa of Percent Moisture by Vol>«ne 

144: on, )* 0*0(11 j0!31*R*0(3) ; 0 ( 25 ) • F*0 ( 25) j 01 4 2) *3 *0 ( 42 ) 
145: F’(l)*0*Pni ;P(ai*e*P(Ul ; P1251*F*P(251 jP(121*3*P(421 
*23434 


Calculation of Percent Air 

146: l-om-W[2)*’jm 
147: 1-P(l)-W(21*vm 
143: l-0(3)-W(31*Ul31 
149: 1-P(81-W|31*V(3) 

113: 1-01251-;J(41*:.'(251 
151: 1-P1251-;M4)*V(251 
152: 1-D1421-W15) *U(421 
153: 1-P1421-W15)*V(421 
*17723 


194 • 


151i 

lS5t 

15C( 

157« 

lS3i 

iy*t 

160> 

16It 

162: 

153: 

154: 

155: 

165: 

157: 

150: 

159: 

*79)3 


Of Pouf. .. -■ 

Specification of Thcrfflal CoiiUuctivitic& 
l*I:2*Xjrll 'CC:iauc A' (U 1 1 1, 0 H 1 /•MX) ) 

(0(l)M;lll4.rl*j|2)*K(2)+r2*U|ll*Kl0))/(Otl)-Hl*'.M2}4-r2*U(l|)-C|l) 
l^I;2»X:cll 'CO-.DUC D' (V ( I 1 , P [ I J ,;i| X) ) 

(P[1)*K{1! *K(3)*r2*V{l)*K|0! ) / (P( 1) *r3 *.;) 2) +c2‘ V) I )) |J ) 

8*I;3*X? rll 'CONDUC A' {U ( I ) , 0 f I ) /M X ) ) 

(3|3)*.<U)*rl*VM31*K|2I+r2*Ut31M((01 )/(Of8) + rlMM31+r2*U(Cn *C13J 
8*i:3*xjcli 'coriouc n' {vf 1 1, pf i ) ,;M'<) ) 

(PI3)* <m ♦C3‘‘.M3)*K 13 )+r2*V(«)*K|0) )/(Pl81 + i3*;M3)4r2*V{9) )*D(3 I 

25-I;4*X;cU 'COIIOUC A' (1) 1 1 ) ,0 ( I ) /M X) ) 

<3(251*i;(l)fcl*',M4 )M:t2)-*t 2Mjf25)*K10) )/(0(25)4rl*;M'M+i2*U(25) 

25-1: 4-xjcii 'cotiauc n' (V (1] ,p( 1 1 /M.xn 

(P(25)*Ulllfr3*W(4 )*K(3)-^r2*V(25|*lM0) )/ ( P [ 25) +c3* ( 4 ) * V [ 25 ) ) -0 ( T*; 1 

42*I;5*XicU 'CO'IDUC A' ('J ( I ) ,0( I )’.MX ) ) 

Ot-**2)*Kll)-: clS'I5)*:{(2)+c2^'J(42)*(c;on/(v')(42)4rl*W(51+r2*U(4 2n •Cf12J 
4 2>lr5*Xjell 'co:iouc »J ' (V {1 ) ,P{ I ) , .M X)) 

(P(4 2)*Xlll4-r3*W(51*Kt3)+r2*V(42) *K(0| )/(P(421+r3*Vi(5)+r2*V(42M-D(-i2) 


Specification of Heat Capacity 

170. V}(2)‘C:2)tOm»'^fl]*En I 
i:i: W12) *3l?)4 P(U» .lll*FUl 
172: IM31 "G(21+0(31 Vtm *E(a) 

173: ;M3) •3(2|4-Ptni *:MI i*Ff 31 
174: .M4) 3(21f.->(25»‘3(l) *E|251 
175: W|41 •::f2)*P{251*3|l|4^?»25) 

176: W|3)*G! 2) Oiw; *3 (1! -Cf421 
177; kMS) *C! 2l4.pH2)*a (1) *F142J 
*167x3 


Interpolation of Thermal Conductivities and Heat Capacity 

173: Cot 1-2 to 24 jC(91-{C(25)-C{3n/17M3-I)>ClIj 
179: 0(31-(0(2S1-D{!^| )/L7* ) 

100: i2[2l - (S|251-':lPl )/l7*(3-I)-VJ(I I 

101: Pl3)-(:-'|251-K(31)/17M3-I)*c’fI)jnaxt I 

182: Cor :«26 to 50 ;C f 25 1 (C |42 ) -C ( 25 1 ) /17* ( I -25) -C f I ) 

1S3: 0(25]+{D{421-ni23))/17* (I-25 )*t(I) 

184: Gf25I + {E(421-3(251)/17Ml-25) •Eli) 

18S: F(25j4-{F|42)-E(25n/17* (I-25)*tMI ) >next I 
165: C{l)-Cf01 ;n[ll^Dl01;E(l)-S(01;F(ll-F)0) 

•4021 
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Prjintout of Thermal Conductivity (cal/®C cm) 

lQ7t fnt l,7/,'C.')'JVJCTlVIYy"fwrt 706.1 

108i fmt i,3/,4x,)"OCPTH SITB A SITE D "|wct 706.1 
I09i lnt(:l/3»^l*irk 
190; c2*3-c3 
l«l: M-c3^l*c4 

I92i for 1*3 to c2-l>I*r5*r6 
I93t if r4»l:E64>l>c6 
194: r6*c2*r7*ta 
195: if r4>>22 ; r3f l«c3 
196: r0*r2*r9 

197: fmt l,5x,f3.0,3x ,f7.5,2x,f7.5,4K .f 3.0,3x,f7.5,2x,C7.5,ss 

193: wrt 705 . 1 , rS ( ri J, D( r5 ) , r 7. C ( 1 7 ) , D ( r 71 

199: fnt 2 , 4x , f 3 . 0 , 3x , f 7 . 5 . 2x , f 7 . S 

200: wrt 706. 2« r9 .C ( r9 h 0 [r9 1 ; next 1 

*31900 


Printout of Volumetric Heat Capacity (cal/cm^ sec ®C) 

201: fnt 1, 9/, “HEAT CAPACirr'jvrt 706.1 

202: fmt 1, 3/, 4x , 3 "DEPTH SITE A SITE B "jwrt 706.1 

203: int(H/3) tl«*r2 

204: r2*3*r3 

205: N-r3+l*c4 

206: for l»3 to r2-lfl*t5«r6 
207: if r4>»i ; r3+l^c6 
238: r6-i’r2*c 7*r3 
209: if r4>-l; r6+l*r0 
210. r8+r2*r9 

2)1: fmt l,Sx,f3.0«4x,fC.3,3x,f6.3.4x .f 3.0,4x,f6. 3«3x,f6.3,z 
212: wrt 706. 1, rS , E (rS | , r' ( r 5 1 , r 7, E ( r 7 1 , F( 1 71 
213: fmt 2, 4x , f 3 .0 , 4x , £6 . 3 ,3x , f G . 3 
214: wet 706.2,c9,E|r9),?Ir9);next I 
•16222 


Calculation of Thermal Inertia 

215: for K«0 to 50 ■ 

216: (C(K1*E(K) r.5-I|Kl 
217: {DIK1*F(K1) 

*24481 


Calculation of Thermal Diffusivity 

218: ClKl/GlKl-OCvl 
219: D t«I /F1141*TIK1 
220: next K 
*29332 
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Printout of Thermal Inertia 

221: Cmt I, 9/, " jwrt 706.1 

222: Cmt I, 3/, lx ,2 "DEPTH SITE \ SITE S “>wit 706.1 

223: int(:j/3)+X*c2 

2Z^^ c2*3*c3 

225: N-c3 + l*r'. 

226: for I»0 to r2-l;I*r5*c5 
227: iC r4>-L>rO>l*i6 
22P : r6*r2*c7*r 0 
229: if r4>-2> c3 + 1t 6 
230: r3*r2*r9 

231.1 Cmt l,5x,£3.0,3x.C7.5,2',f7.5,4x ,£3.0,3x,C7.5,2x,C7.5,z 
232: wrt 706.1,r5,I|r5l,J!r51,r7,l|r7],J(c7) 

233: fnt 2, 4x , f 3.0 , 3x , f 7 . 5, 2x , C7 . 5 
234 : wi t 706. 2,r9 , 1 (r9 ),J (r9 ) ;ncxt I 
•13408 


Printout of Thermal Diffusivlty 

235: fmi- 9/, "THCR'UL. UIFFUH9 I VI7Y" ; wr t 706.1 

236: Cnt ,V . 3 /, ; x , 3 " DEPTH HITE A SIPF. P ";wrt 706.1 

237: int( l/3)+l»r2 

233: r2*3^r3 

239 : N~c3H*r4 

240 : for I*»0 to r2-l j I-*r5*r6 
241: ir r4>»l! c0 + 1t6 

242 : r6+r2-*r 7-*i'3 
24 3: i £ r 4>*2 ; i3 + l*r9 
24 4 : c3-*r2»r9 

245: Imt 1 , 5x , C 3 . 0 , 3x . f 7 . 5 , 2x , C 7 . 5, 4x , C 3. 0 , 3x .£7 . 5 , 2x , £7 . 5 ,z 
246 : wrt 70G . 1, r 5 ,;>! r 5 1 , T ( r5 ) , r 7,0( i 7 ), T ( c7 ) 

24’: Cmt 2,4x,£1.0,3x,f7.5.2x,£7.5 
2i3: wrt 706.2, t9,0Ir91,T(r3)jncxt 1 
•25223 


Time Interval Between Calculations 

249: T/(2»X(ll*XUn*Cl jgio 2f>4 
•551G 


Call for So.il Heat Flux Subroutine 

250: cll 'OHEAT' 

*11251 
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Calculation of Nodal Temperatures and Lower Boundary 

25lt tot I-l to !J-1| (njn/rI-C(l-l|-2*C|n-C|Uin*Am*rO 

252 1 ((cii-n^cm)*Aii-iiMcm+cimn*A!i+i)*-c6)*ri/Gm*nm 

253 » (Fm/ti-on-ii-2*Dm-o(i*n )*»in*t7 

254* (' 0|I-ll^0( II )* 311-lJ* (Oil )^0|I + l))*n(X + U+t7) *rl/Fin •cm ;next X 
255: 4*xm*tlU(Em/cl-2*Cm-2*CUn*M0l4.(2*CUl+2*C|2M*MH*r;(01 
25G: rl/nm*niO) •niOl 

257: 4*xm*rl2M Fill /r 1-2*0111 -2*1) l2n*O|0l4(2'O|in2*D| 21) *om •SI 0) 
25St rl/F|ll*S(31*S|01 
250: Mi]l«n(.M 

260: Q(:n«s[:n 

•32346 


Reassignment of Nodal Temperatures for Succeeding Intcration 

2C'i! foe l«0 to :4:R(I1*MI1 rSIH^aUltnoxt I 

262: M01>n(01«'MSll 

*29200 


Test for Printout Time 

2b3; if Ml41<l‘jgto 305 
•256 


Printout of Pertinent Data 


264: 

rl0*-l*tl0 


265: 

int(rin/2)*rll jrll*2*r) 2 


2G6: 

fnt 1, 4/,"T')PCn\PURr: ’1DFILB AT ",£z2.0,f7.2 

.0," HOURS 

267: 

wrt 70b.l,.im ,:5(11 


269: 

fnt 2,/,/,"R:)R^^CS TeMPERNTURf. OIFFE»R;JCF, •" 

,f6,3 

269: 

wrt 706.2,3151) 


270: 

fnt 4,/, "SOIL H':A7 flux SITF A • ",£12.9," 

5ITF. 0 

271: 

wet 706. 4, A, 0 


272: 

£mt l,2/,4x,3"OEPTH TEMP A TEMP 3 " 


273: 

wet 706.1 


274: 

lnt(;>/3)+lT2 


275: 

r2*3*r3 


276: 

JI-c3<H*e1 

. 

2'»7: 

foe I«0 to t2-l»I^t5*e6 


: *6; 

if r4>«l : r5*l*c5 


279: 

cS+e2-*r7*c3 


230: 

if e4>«2 : 


291: 

t0+r2*r7 


282: 

fmt l,Sx.£3.3,4x,f6.3,3x,f6.3,4x ,£3.0,4x,f6. 

3,3x,C6.3, 

233: 

wet 705.1,c5,Mrjl,3fr51,c7,Mc7) ,3lc71 


284: 

fnt 2,4x,f?.3,4.{,fr*.3,3x,fr,.3 


285: 

we t 706. 2, c9 , AlrO 1 ,31 c9 1 inext I 


•6549 



fl2.9 


Procedural Step (Reset) 


236: 0*'1(41 
•3S74 
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Call for ^’lot Routine 

2S7t if rl2lrlfl;qto 3C6 
2391 cU 'Plots' 

•19455 


289 

2:>o 

2n 

292 

293 
29‘. 
293 

296 

297 
293 


Plot of Theoretical and lixpej ir.cr.’tal Soil Tcrijn'ratures 

for I"0 to 'iD»on 'olotA' (All) ,L(Il)>ncxt I 
if !|(H-3)24»‘!(l I 

for J’l to Sjcll 'roalA' (ZIJ ,J!|11/2,J.) ,l}(3n ;next J 

for I "O to 50;c)l ' olotJ ' |fi 1 1 I , L ( I 1 J j Ncx f i 

Coe .T-1 to j;cU 'cei]A' (2(J,!l(l |/:,2’ ,:J(.n ) ;ncxt .7 

fnt “:wct 706.1 

fot 4,/,"3rrr. A-,zjwct 700.4 

f rt 2; 3x,C3.0,3;< , C3.2, 4x,C3.0, 3< ,f 5. 2,4x ,f3 .0,3x »E3. 2,z 
fmt 3.4:;.f3.0,Jx,f3.2,4x,f3.0,3x ,f S.2 
wrt 703 


Plot Of Theoretical and Experimental Soil Temperatures 

239 5 Cor I-'O to SOjcll ' olotA' (M I ) , L {I ) ) ; next I 

290 5 if 'l|U« 0 j /.4*:tfl) 

2915 Co.- d»l to 5jcU 'roftlA' (.: 1 7 ,?M 7) ) jnext .7 

2S2: for I-O to SOscll 'olot3 ' (3 , 1 ) .L ( I ) ) 5 nox t I 

293: for J»1 to Sscil ' t calB ' (2 [J , H ! 1) /2, 2) l.l 1 > 5 next J 

•6936 


Printout cf Experimental Soil ‘^^'mperatures for Piofae A 


29^: r.nt 
295s fr.,;: 
796: f;'t 
297: Crd 
;'9as vrt 
2 >9 5 wrt 
•2683 


1 , / ,9r ,5''L»r.rni TCMP ";wrt 706.1 
4, rs A'* ,25 wrt 706.4 

2,3x,i3.0i .'^t,C5.2|4x,f3.0,3X|C5.2,4x,. i.0,3''tr7 j.2,7 
3,4\,C3.0,3x ,f5.2, 4x,E2.0,3x ,C5.2 

70 >. 2 ,.:[l) , 2 ( 1 , 3 ( 11 / 2 . 1 ) ,:J( 2 ) , 3 ( 2 , Ml n/ 2 , 1 ) ,:J 1:1 , 1', ! '! l! ( 1 1 . i! 
706.3,;U41 ,2[4,um/2,l) , 1 J ( 5) , 2 I 3 , 11 { 1) /2 , J. ) 


Priatovt of E.\porimental Soil Temperatures for Profile; 1\ 

300s Cmt 5,/,"sn’E Cf.aswr'; 706.5 

3C3 5 fT.t 2,3x,f 3.0,3x ,C5.2, 4». ,f. 3.0,3x ,f 5.'/’ , 4x ,f 3.0,3x ,fS. ?,s 
:;0:s fr.t 3,4x,C3.0,;x,C5.2,4/,£3.0,3x ,C5.n 

333 s wrt 706. 2,:i(H ,2(1 ,H (11/2, ?.1 ,:i 121 , Z C,M (11/2,2) ,;J (3! , 2(3,!! I. M/2,2) 

304 t vrt 706. 3,N (5 1 , 2 [ 4 , 3| D /2, 2) ,‘l ( 5 ) , 2 (5, H ( 1.) /2, 2 1 

•24421 


Test if Calculation Has Run Desired Time 
30S. if H(l)-24 ;0-*ll(H 

306: ii MUI<ri(21 sntn 313 

307: if H(lj<;!f Jl ;qto jlO 

309: if Yl0J<y (1) solo 310 

30 9: oto 32 4 
‘711 


Ttst If Tmpomure Difforonct it to bo Ctlculated 

)iOi 

)Ui 1C M3|/(cl3*20)>«li9to 313 
3I2i qto 3L7 
3l3t rl3«l*cl3 
•30522 


Calculation of Surface and S cm Tenperature Oiffcrenco 

314: M0)-3t01*-|!1I3I/20l 
3ISi M5)-D(S)«Vt;i(3)/20i 
•30792 


Calculation of Navr Time 

3I6t M(31/60*0(M(3|/30) 

3l7i nl4I^r*Ml41 
^I8i MUJ^T/eO.Mm 
319t it M(ll< 60 ) 9 to 323 
320t 

321( :iC Hm<24;qtO 323 

322: liUt-24«tt(\) }tt0Ul«Y{0) 

323: 9 to 14 L 
324: oto 
•28526 


Plotting of Surface and 5 cm Temperacuro Differences 
325: ell 'OifC oLot' 

326: for I«L to 72:ell 'olotS' (0(11 «UtXI) :ncKt I ’ 

327: tor I*l to 72jcll ' plotC ' |0( 1 1 , V( 1 1) i next I 

320: end 

•29769 


